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SUMMARY 

Laser  propagation  through  the  atmosphere  is  severely  limited  by 
atmospheric  effects  such  as  fog  and  rain.  Because  of  the  high  power 
achievable  and  convertness,  CO^  lasers  are  finding  a multitude  of  uses 
in  military  applications  such  as  target  designators,  be?”'  riders,  moving 
target  indicftors,  terrain  following,  weapons  delivery,  and  imaging. 

CO^  lasers  have  also  grown  in  importance  in  applications  such  as  pollution 
monitoring,  optical  communications,  and  earth  resources  remote  sensing. 

These  applications  require  a detailed  knowledge  of  CO^  laser  atmospheric 
extinction  and  its  dependency  upon  relevant  meteorological  parameters 
associated  with  adverse  weather  phenomena. 

A tunable  CO^  laser  transmissometer  with  a 12  dB  dynamic  range  was 
fabricated  to  measure  the  extinction  of  laser  radiation  under  adverse 
weather  conditions.  Data  on  eight  different  lines  has  been  collected 
using  a phase-lock  receiver  with  an  optical  aperture  three  times  the 
beam  diameter.  The  system  dynamic  range  of  72  dB  is  limited  by  the  use 
of  pyroelectric  detectors.  It  was  found  that  during  dense  tog  measure- 
ments, the  extinction  encountered  exceeded  the  dynamic  range  of  the 
system  more  than  50“  of  the  time.. 

Reported  measurements  include  extinction  as  a function  of  wavelength, 
time,  and  atmospheric  particulates.  Experimental  measurements  in  fogs 
with  extinction  ranging  from  1 to  72  dB/km,  in  rain  with  extinction 
ranging  from  2 to  17  dB/km  and  in  light  snow  wi '■.h  extinction  ranging 
from  3 to  14  dB/km  are  reported.  Partic^p  size  data  for  over  150  fog 
size  distributions  have  been  collected.  Extinction  by  fog  in  excess  of 
290  dB/km  was  obtained  from  computed  values  of  extinction  cross-sections 
and  measured  size  distributions.  An  empirical  m'^del  was  derived  for  CO^ 
laser  extinction  as  a function  of  precipitation  rate.  Depending  on  the 
refractive  index,  e.itinction  by  haze  particulates  may  exceed  1.4  dB/km. 
Particles  below  0.3  um  radius  do  not  contribute  significantly  to  extinction, 
and  unless  the  aerosol  size  distribution  contains  a significant  number  of 
particles  larger  than  this  size,  aerosol  extinction  may  be  neglected. 
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Measurements  show  that  transmission  losses  in  dense  fog  can  easily 
be  six  orders  of  magnitude  greater  than  losses  in  heavy  rain.  Results 
indicate  that  extinction  of  CO2  laser  radiation  in  adverse  weather  is  at 
best  a very  weak  function  of  wavelength. 
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SECTION  I 
INTRODUCTION 

Laser  propagation  through  the  atmosphere  is  severely  limited  by 
atmosphere  effects  such  as  fog  and  rain.  Because  of  the  high  power 
achievable  and  convertness,  CO^  lasers  are  finding  a multitude  of  uses  in 
military  applications  such  as  target  designators,  beam  riders,  moving 
target  indicators,  terrain  following,  weapons  delivery,  and  imaging.  CO^ 
lasers  have  also  grown  in  importance  in  applications  such  as  pollution 
monitoring,  optical  communications,  and  earth  resources  remote  sensing. 
These  applications  require  a detailed  knowledge  of  CO^  laser  atmospheric 
extinction  and  its  dependency  upon  relevant  meteorological  parameters 
associated  with  adverse  weather  phenomena. 

The  object  of  this  study  is  to  investigate  how  atmospheric  particu- 
lates, especially  fog  and  rain,  affect  CO^  laser  propagation  through  the 
atmosphere. 

Atmospheric  particulates  below  0.3  mictometer  radius  add  little 
contribution  to  extinction  at  CO^  wavelengths  and  generally  can  be  ignored 
in  computations.  In  dealing  with  hazy  conditions,  most  of  the  aerosol 
contribution  to  extinction  comes  from  particles  in  the  size  range  0.3  to 
1.5  n”!  radius,  with  smaller  contributions  extending  to  particles  of  radius 
up  to  4 pm.  For  typical  summer  weather  under  hazy  conditions,  molecular 
extinction  due  to  the  far  wings  of  the  water  vapor  lines  predominate, 
with  an  approximate  upper  bound  contribution  to  extinction  of  .34  km"^ 
due  to  aerosols  (this  locality,  Dayton,  Ohio). 

The  fog  aasurements  presented  are  backed-up  by  an  f.'u-osive  particle 
size  data  base.  Computation  of  extinction  by  fog  has  beer,  performed  by 
direct  laser  measurements  and  by  numerical  integration  of  the  extiOv+ion 
cross  section  and  particle  size  distributions.  The  liquid  water  content 
for  each  di'-tribution  has  been  computed  and  this  is  used  to  obtain  an 
approximation  to  the  extinction  which  is  then  compared  with  the  measure- 
ments. Extinction  due  to  rain  is  presented  as  a function  of  rain  rates 
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and  a mathematical  relationship  is  derived  relating  extinction  to  rain 
rates  without  the  need  for  knowledge  of  rain  drop  size  or  drop  terminal 
velocity. 

The  principles  of  molecular  spectra  are  presented  in  Section  II, 

The  derivations  and  formalism  covered  are  centered  around  the  two  main 
contributors  to  molecular  extinction  near  the  10.6  ^!m  - the  H^O  and  CO^ 
mol ecul es . 

Molecular  extinction  is  treated  in  Section  III.  Line  intensity 
profiles  and  wing  theory  are  treated,  and  an  introduction  to  high  reso- 
lution transmission  models  is  presented. 

Section  IV  discusses  extinction  by  atmsopheric  particulates  and 
precipitation.  Various  model  distributions  are  presented.  Exact  Mie 
extinction  formulas  are  introduced  as  well  as  the  Deirmendjian  approxi- 
mations . 

Atmospheric  turbulence  effects  and  how  they  affect  system  design 
are  presented  in  Section  V.  Analysis  of  the  laser  transmissometer  system 
is  presented. 

Section  VI  describes  measurements  and  analyzes  the  results  obtained 
for  CO^  laser  extinction  in  fog,  rain,  and  haze. 

The  data  collected  for  this  study  stretches  over  a three  year  period.. 
This  study  depended  on  the  ability  to  obtain  the  proper  weather  conditions 
of  interest.  Sometimes  when  the  weather  was  right,  the  equipment  was  not 
ready;  and  when  the  equipment  was  ready,  the  weather  was  not  right. 

Several  fog  measurements  were  obtained  during  the  winter  of  1979,  but 
the  particle  counter  data  was  suspect.  Fall  and  winter  of  1980  produced 
only  one  brief  but  unexpected  fog.  The  results  presented  in  this  report 
are  based  on  data  collected  within  the  last  year,  mostly  because  the 
year  1981  provided  all  the  necessa-"'  adverse  weather  conditions  of 
interest,  and  calibration  of  the  meteorological  instruments  needed, 
achieved  a good  level  of  confidence  which  was  absent  in  previous  years. 
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Several  difficulties  were  encountered  during  the  course  of  this  in- 
vestigation (See  Appendix  C)-  Troubles  with  the  laser  were  many.  The 
two  most  persistent  ones  were  the  line  selection  mechanism  and  gas  leakage 
from  the  plasma  tube.  The  mechanical  design  for  the  grating  drive  made 
it  impossible  to  maintain  calibration  on  the  laser  for  more  than  a few 
days.  Gas  leakage  of  the  plasma  tube  resulted  in  the  laser  exploding 
on  three  different  occasions.  Because  of  the  unreliability  o^  the  large 
tunable  laser,  a smaller  single  line  laser  was  modified  and  fitted  with 
a piezoelectric  mounted  mirror.  This  arrangement  provided  data  on 
approximately  ten  lines,  all  of  which  were  not  available  at  all  times.; 
Where  data  on  specific  lines  is  absent,  it  is  because  the  line  could  not 
be  tuned  in,  or  because  it  took  longer  than  required  for  under  conditions 
of  the  measurement. 

In  this  report,  extinction  is  defined  as  the  loss  of  radiant  flux 
(or  power)  along  the  direction  of  propagation  due  to  the  interaction  of 
electromagnetic  radiation  and  matter.  As  defined,  it  is  the  combined 
contribution  of  molecular  and  aerosol  extinction.  No  effort  is  made  to 
separate  scattering  effects  from  absorption  effects.  The  extinction 
coefficient  is  the  value  of  extinction  as  measured  or  computed  at  specific 
frequencies  (or  wavelengths)  per  unit  path  length. 
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SECTION  II 
MOLECULAR  SPECTRA 

The  extinction  of  an  electromagnetic  beam  by  an  absorbing  medium  is 
a result  of  the  inte'^action  of  the  electromagnetic  fields  of  the  radiation 
with  the  molecular  species  in  the  absorbing  medium.  If  the  interaction 
results  in  the  molecule  taking  energy  from  the  field,  absorption  occurs. 

If  the  interaction  results  in  the  molecule  giving  energy  to  the  field, 
emission  occurs.  The  molecular  species  emit  or  absorb  at  specific  fre- 
quencies which  depend  on  v^eir  molecular  structure.  At  the  frequency, 
where  a molecule  emits  or  absorbs,  it  is  said  to  have  an  absorption  line. 
Thus  an  absorption  (or  emission)  process  results  in  an  absorption  spectrum 
composed  of  discrete  lines,  each  corresponding  to  an  energy  state  as  pre- 
scribed by  quantum  mechanical  principles.  The  absorption  of  energy  occurs 
when  the  atomic  system  makes  a transition  between  one  stationary  state  of 
energy  E'  and  another  of  energy  E",  where  by  convention  the  lower  of  the 
two  levels  involved  in  a transition  is  indicated  by  a (")  and  the  upper 
of  the  two  levels  will  be  indicated  by  a (').  The  frequency  at  which  the 
transition  occurs  is  given  by 


h 


where  h is  Planck's  constant. 

The  appearance  of  a spectrum  is  connected  with  energy  alterations 
in  the  molecular  species  which  produce  the  spectrum.  Molecules  can  uptake 
energy  in  four  possible  ways: 

1.  Electronic  energy  (E^) 

2.  Vibrational  energy  (E^) 

3.  Rotational  energy  (E^) 

4.  Translational  energy  (L^). 
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The  translational  energy  is  riot  quantized  and  does  not  result  in  or 
influence  spectra.  The  total  energy  which  produces  the  spectrum  is  given 
by 


E 


= E. 


a 


+ Ey  + Er  + El 


(2) 


where  E^.  is  the  interaction  energy  between  the  different  energies  in  the 
molecule.  Each  of  the  energy  states  is  described  by  certain  quantum 
numbers,  an  alternation  of  which  leads  to  a specific  spectrum.  As  these 
spectra  are  found  in  different  regions,  the  energy  changes  connected  with 
the  emission  or  absorption  of  radiation  can  be  of  different  orders  of 
magnitude.  The  largest  energy  values  are  found  for  the  energies  of  electron 
promotion.  The  allowed  electronic  energy  states  are  rather  widely  separated 
in  energy,  so  that  transitions  between  different  energy  levels  correspond 
to  transition  frequencies  which  lie  generally  in  the  visible  and  ultra- 
violet portions  of  the  spectra.  These  electronic  bands  play  a very  minor 
role  in  the  transfer  of  radiation  through  the  atmosphere  and  they  are  of 
no  consequence  when  dealing  with  absorption  spectra  associated  with  CO^ 
infrared  lasers. 


Interpretation  of  the  spectrum  of  a molecule  requires  a knowledge  of 
the  allowed  energy  states,  allowed  transitions,  and  of  the  transition 
probabilities,  and  how  each  of  these  depend  on  the  molecular  structure. 

To  understand  infrared  spectra,  one  must  know  how  vibrational  and  rota- 
tional motions  of  a molecule  produce  their  associated  energies.;  A complete 
treatment  of  molecular  vibrations  and  rotations  requires  the  application 
of  quantum  mechanics  and  group  theory.  These  methods  become  cumbersome 
and  impractical  even  for  simple  molecules.  Knowledge  of  the  molecular 
structure  makes  it  possible  to  apply  symmetry  properties  which  allow  seme 
simplification  in  arming  at  a solution.  In  general,  the  most  common 
approach  requires  the  separation  of  the  rotational  effects  from  the 
vibrational  effects  and  neglecting  the  cross-coupling  effects.  Higher 
order  approximations  necessitate  the  inclusion  of  vibration  interactions 
and  rotation-vibration  interactions,  where  the  mathematical  formulation 
has  been  made  to  agree  with  observed  results. 
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A vibrating  or  rotating  molecule  cannot  interact  with  electromagnetic 
radiation  unless  the  vibration  or  rotation  is  accompanied  by  a dipole 
moment.  An  oscillating  dipole  must  be  present  for  coupling  to  occur 
between  the  energy  of  the  molecule  and  the  electric  field  so  that  energy 
can  be  exchanged  between  them.  This  restriction  precludes  all  homonuclear 
diatomic  molecules  such  as  0^,  and  N^,  from  possessing  an  infrared 
spectrum.  Under  normal  conditions,  their  spectra  can  only  be  obtained  in 
connection  with  electronic  promotion.  Therefore,  in  the  study  of  CO^ 
laser  propagation,  it  is  found  that  CO^  and  H^O  are  the  primary  molecular 
absorbers  in  the  lower  atmosphere..  In  the  discussion  that  follows,  the 
vibrational  and  rotational  energies  of  these  two  types  of  molecules  will 
be  derived,  showing  how  they  give  rise  to  their  absorption  spectra. 


For  a molecule  containing  n atoms,  the  one-dimensional  Schroedinger 
equation  can  be  easily  extended  to  describe  the  system  under  consideration. 
For  a system  of  N particles,  the  wave  equation  is  given  b' 


_1 + Vi|i  * E\p 


(3) 


where  m.  is  the  i-th  particle  mass,  E is  the  energy  operator  (-h/Ziri) 

2' . 

3/3t,v.j  is  the  Laplacian  for  the  i-th  particle,  and  the  wave  function 

,y,  ,z, , . . .X  ,y  ,z  ) involves  all  3N  cartesian  coordinates.  In 
111  n n n 

general,  three  cartesian  coordinates  are  required  to  describe  the  motion 
of  the  center  of  mass  of  a molecule.  These  are  used  to  describe  the 
vibrational  modes  of  the  molecule.  Three  more  coordinates  are  required 
to  described  the  rotation  of  the  system  about  the  principal  moments  of 
inertia  axis.  Wilson,  Deciuss,  and  Cross  (Reference  1)  have  shown  that 
the  system  can  be  treated  separately  for  each  set  of  coordinates,  so  that 
the  vibrational  motion  can  be  separated  from  the  rotational  motion.  If 
the  interaction  energy  terms  are  neglected,  the  total  wave  function  can 
be  written  as  the  product  of  a vibrational  wave  function  and  a 
rotational  wave  function  Yp,  i.e,. 


6 


AFWAL-TR-81-1280 


A molecular  system  composed  of  N atoms  having  3N  coordinates 
possesses  3N  degrees  of  freedom.  Those  degrees  of  freedom  include  three 
that  describe  translational  motion  of  the  system  as  a whole.  If  one 
chooses  the  molecular  center  of  gravity  as  the  origin  of  the  coordinate 
system,  then  the  total  number  of  coordinates  is  reduced  by  three,  since 
3N-3  coordinates  are  sufficient  to  determine  the  relative  position  of 
the  atoms  with  respect  to  their  center  of  gravity.  Additionally,  since 
the  vibrational  and  rotational  wave  functions  are  separable,  the  number 
of  coordinates  is  reduced  by  3 for  either  a rotational  or  vibrational 
system.  Therefore  3N-6  coordinates  are  needed  for  a description  of  a 
purely  vibrational  system  or  a purely  rotational  system.  Due  to  symmetry 
considerations,  a linear  molecule  has  only  two  degrees  of  rotational 
freedom;  therefore,  3N-5  degrees  of  freedom  are  required  to  describe  a 
linear  molecular  system. 


1,  MOLECULAR  VIBRATIONS 


The  quantum  mechanical  calculation  of  the  possible  vibrational  fre- 
quencies for  a given  molecular  system  starts  with  the  solution  to  the 
wave  Equation  3.  In  cartesian  coordinates,  the  wave  equation  is 


i=l 


1 

( — + + 


) + V(x.,y 

1 i 


• • 2^)  i|) 


(4) 


• Eij; 


where  the  position  of  each  atom  of  mass  m.  is  described  by  the  triplet 
(x^.  ,y^  ,z^. ) , The  use  of  cartesian  coordinates  is  not  convenient  for 
describing  molecular  motions  and  one  resorts  to  the  use  of  normal  coor- 
dinates Q|^,  k = 1,2,... .3N,  where  the  normal  coordinates  are  defined  in 
terms  of  the  mass-weighted  cartesian  coordinates  q^  by  the  linear  trans- 
formation 

N 

\ “ Zj  ^ * 1.2,...3N  . 

i=l 
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Normal  coordinates  are  such  that  both  the  kinetic  and  potential  energies 
of  a system  are  expressed  by  the  sum  of  terms  each  involving  the  square 
of  a normal  coordinate.  Pauling  and  Wilson  (Reference  2)  have  shown  that 
in  terms  of  these  normal  coordinates  the  kinetic  and  potential  energies 
are  given  by 

3N 

T - % IT 

k-1 


V 


k-1 


2 

k'^  k 


XvQ 


(7) 


where 


(8) 


represents  a change  of  scale  of  the  general  coordinates 
designates  the  position  of  k-th  nucleus  with  respect  to  its  equilibrium 
position,  is  the  mass  of  k-th  nucleus,  and  are  the  constants  related 
to  the  normal  frequencies  X.  In  terms  of  the  general  cartesian  coordinates 
the  first  term  in  Equation  4 can  now  be  expressed  by 


By  the  change  of  coordinates  of  Equation  8 

_1l_  _J_J! — becomes  — 

«k  si'k  ’'?’k 
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and  after  employing  Equation  3,  the  wave  equation  (Equation  4)  can  be 
written  in  the  form 


(9) 


Using  the  fact  that  for  a purely  vibrational  system,  the  system  can  be 
entirely  described  by  3N-6  coordinates  {3N-5  for  a linear  system  of  atoms) 
and  rearranging  Equation  9,  the  vibrational  wave  equation  becomes 


3N-6 

3N-6 

S 

k=l 

3 + 8^-1  E - V 

E)  - 0 

k*l 

(10) 

• «i.  1 ^ ^ 

3Q^k  ^ L 

Letting 

<1^  “ 'l'i(Ql)'i'2(Q2^ 

•••'^3N-6^^3N-6^ 

(11) 

Equation  10  is  now  separable  into  3N-6  equations,  one  for  each  normal 
coordinate  Qj^i  each  one  having  the  Form 

= + (12) 
k ^ 

This  is  the  one-dimensional  equation  for  the  simple  harmonic  oscillator.. 
The  solution  (j;  of  the  vibrational  problem  is  therefore  expressible  as  a 
product  of  3N-6  harmonic  oscillator  functions  while  the  vibrational 

energy  E is  the  sum  of  the  energies  of  3N-6  harmonic  oscillators,  and  is 
thus  given  by 

E . + Ej  + ...  -f  (13) 

The  energy  levels  of  a linear  harmonic  oscillator  are  given  by  the  well 
known  relation 

E = (v  + %)hv  V “ 0,1,2,...  (14) 
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where  v is  the  quantum  number  which  can  take  on  any  positive  integral 
value,  while  v is  the  classical  frequency  of  the  vibrating  system  and  is 
related  to  the  wavelength  by  the  relation 


X 

V * 

2ir 


(15) 


Using  Equation  14,  the  vibrational  energy  of  a molecule  with  several 
classical  frequencies  corresponding  to  the  normal  modes  of  vibration 
is  given  by 

M 

E(v)  » E (v^  + %)hv^  (16) 

i“l 


where  M is  equal  to  3N-6  or  3N-5  for  the  case  of  linear  molecules.  When 
two  or  more  normal  vibrations  are  equal.  Equation  16  has  to  be  modified 
to  account  for  the  degree  of  degeneracy  (d^. ) and  the  result  is  given  by 

M 

E(v)  “ ^ (v^  + d^/2)hv^  (17) 

i*l 

Where  degenerate  terms  occur,  the  energy  value  will  depend  not  on  the 
individual  values  of  v^  and  V|^  but  on  their  sum.  If  both  sides  of 
Equation  17  are  divided  by  he,  then  the  energy  term  is  given 

M d 

G(v)  - ^ 

i-1 

where  both  6(v)  and  u.  have  the  wave  number  units  of  cm’\ 

Equation  18  gives  a first  order  approximation  to  pure  vibrational 
spectra.  In  the  study  if  pure  vibrational  spectra  of  polyatomic  molecules, 
it  is  found  that  in  addition  to  the  bands  corresponding  to  the  normal 
vibrations,  there  are  bands  of  lower  intensity  present  at  approximate 
multiples  of  the  normal  frequencies.  These  overtones  are  the  result  of 
the  anharmonicity  of  tie  vibrations. 
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When  the  classical  harmonic  oscillator  is  considered,  the  potential 
energy  function  is  treated  as  having  the  form 


(19) 


Examination  of  the  potential  function  for  a molecular  system  shows  that 
the  potential  energy  curve  is  better  approximated  by  the  inclusion  of 
cubic  and  quartic  terms  so  that  Equation  19  takes  the  form 

V - gr^  - hr^  + Rr^  - ...  (20) 


The  consequence  is  that  the  vibrations  are  no  longer  harmonic  in  nature 
and  the  concept,  of  normal  vibrations  is  no  longer  rigorously  possible. 
The  influence  of  this  anharmonicity  shows  u'i  in  the  energy  term  as  an 
interaction  containing  two  or  more  quantum  numbers  in  addition  to  terms 
containing  single  quantum  numbers.  Using  perturbation  methods,  Nielsen 
(reference  3)  arrived  at  the  solution 

3 3 d 

0)^  (v^  + — j^-)  4-  ^ 

i kil 

it 

i k>l 


' T. 

i 


for  molecules  of  the  type  XY^  having  a two-fold  degeneracy,,  "he  last 
quantity  in  the  expression  accounts  for  the  fact  that  some  energy  levels 
are  split  by  the  action  of  anharmonicity,  K being  a new  quantum  number 
having  the  values  1.  = v.,  v.  _ v.  . and  g.,  is  some  constant  depending 
on  the  geometry  of  the  molecule.  The  vibration  spectrum  for  the  CO^ 
linear  molecule  can  be  obtained  from  Equation  21,  For  nondegenerate 
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molecules  such  as  H^O,  1^.  = g^. = 0 and  the  vibrational  spectra  may  be 
obtained  from 

GCvpVj.Vj)  - 21  S S *ik 

1-1  1-1  k-l 

(V,  + -i-)  (V.  + J-)  (22) 

1 2 2 

The  values  of  x.j^  and  have  been  computed  for  the  CO^  molecule  by 

Darling  and  Davidson  (Reference  4).  Similarly,  values  of  x..  for  the 

1 K 

H2O  molecule  have  been  computed  by  Dennison  (Reference  5). 

2.  MOLECULAR  ROTATION 

Solutions  to  the  molecular  rotation  problem  can  be  obtained  by 
setting  the  Schrdedinger  wave  equation,  in  three  dimensions  and  using 
polar  coordinates.  The  straightforward  solution  b)comes  very  complicated 
and  to  simplify  the  solution,  one  has  to  consider  special  cases  using 
symmetry  properties  to  the  molecule.  This  results  in  solving  the  wave 
equation  individually  for  classes  of  molecules  that  fit  into  certain 

symmetry  groups  and  several  solutions  to  the  rotational  wave  equations  i .. 

are  necessary.  Realizing  that  in  absorption  spectra,  the  expression 
for  the  rotational  energy  states  allowed  are  of  prime  importance,  it  is 
much  simpler  to  lesort  to  classical  mechanics  to  obtain  the  desired  energy 
states  .- 

The  rotation  spectrum  of  a molecule  is  mainly  determined  by  the 
magnitude  of  its  moment  of  inertia.  The  atoms  in  the  molecule  are  treated 
as  point  masses,  m^. , so  that  the  moment  of  inertia  of  the  molecule  (assuming 
a rigid  system)  is 

N 

I - ^ . (23) 

i“l 
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where  r^.  is  the  perpendicular  distance  of  the  i-th  mass  f*’om  the  rotational 
axis.  According  to  classical  mechanics,  the  inertial  behavior  cf  such  a 
body  is  completely  defined  by  the  moments  of  inertia  about  three  mutual  y 
perpendicular  axes.  The  corresponding  moments  of  inertia  are  callev.  t‘ie 
principal  moments  of  inertia.  To  determine  the  main  inertia  axes,  ont 
must  have  the  knowledge  of  the  symmetry  properties  of  the  molecule.  Axes 
of  symmetry  are  always  perpendicular  to  main  inertia  axes  and  planes  of 
symmetry  always  perpendicular  to  main  inertia  axes.  When  dealing  with 
atmospheric  absorption  at  CO^  wavelengths,  only  two  types  of  molecular 
symmetry  need  be  considered.  One  is  the  linear  molecule  (CO2)  and  the 
other  is  the  asymmetric  top  molecule  (H^O). 

a.  Linear  Molecule 

In  a linear  molecule  such  as  CO^  the  molecular  axis  is  the  line 
joining  the  nuclei  of  the  individual  atoms.  The  moment  of  inertia  about 
this  molecular  axis  is  in  almost  all  cases  neglected,  being  zero  or  very 
small  in  comparison  with  the  other  two  equally  large  ones.  In  this  case 
the  model  of  the  rigid  rotator  can  be  used  together  with  a correction 
term  for  nonrigidity  where  applicable.  The  nonrigidity  term  is  used  to 
account  for  possible  centrifugal  distortion  produced  by  rotation.  The 
rigid  rotato*  is  illustrated  in  Figure  1,  where  the  center  of  gravity  is 
labeled  CG,  m^  and  are  masses  of  rotator,  and  r^  and  r^  the  distances 
from  the  rotational  axis.  From  Equation  23, 

I « ■*' 

which  reduces  to  the  more  common  form 


where  u 


+ m^ 


is  the  reduced  mass. 
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Figure  1 . Rigid  Rotator 


The  kinetic  energy  of  a system  of  moving  particles  is  given  by 


E - % X 


(26) 


In  terms  of  angular  velocity  w,  the  energy  for  the  system  of  rotating 
particles  is  given  by 


E - % 


2 N . 

% » ^ “1^  i 


where  the  relationship  v = wr  has  been  used.  The  angular  momentum  of  a 
rotating  particle  is  given  by  mvr.  For  a system  of  particles  the  angular 
momentum  is  defined  by 


’’  ■ Z! 

i-1 
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which  is  easily  regrouped  into 

N ^ 

*■  - Xl  (28) 

i*l  ^ 

Bohr's  theory  requires  that  the  quantized  angular  momentum  be  of  the  form 

P - J — . J - 0,1.2,...  (29) 

2r 

where  J is  the  standard  rotational  quantum  number.  It  has  been  shown  in 
various  textbooks  (References  1,2)  that  in  order  to  obtain  results  that 
agree  with  the  solution  of  the  results  obtained  from  the  Schrbedinger 
wave  equation,  it  is  necessary  to  replace  J by  [J{J+1)]  ' in  Equation  29 
In  this  case  the  correct  quantized  angular  momentum  is  given  by 

P - ^(J+I)h/2Tt,  J - 0,1,2,...  (30) 

The  allowed  rotational  energies  are  obtained  by  substituting  Equation  30 
into  Equation  27  and  rearranging  to  give 

+ l)(h/2iT)2  (31) 

^ 21 

From  manipulation  of  Equation  1,  it  is  seen  that  when  energy  is  divided 
by  he,  the  resulting  quantity  has  units  of  wave  numbers  (cm'^).  This 
gives  rise  to  the  spectroscopic  energy  term  F(J)  expressed  in  the  form 

F(J)  - Ej^/hc 

“ (h  /Stt^cI)  J(J+1) 

- BJ(J+1)  cxn"^  (32) 

where 

B » h /Stt^cI  cm'^  (33) 
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is  called  the  rotational  constant.  For  absorption  of  radiation  by  a 
rotational  energy  transition  to  occur,  the  molecule  must  have  a permanent 
dipole  moment  and  there  must  be  compliance  with  the  selection  rule  AJ  = +1 
(aO  = -1  must  be  satisfied  for  emission).  The  wave  number  of  the  radiation 
absorbed  by  an  alteration  in  the  rotational  quantum  number  is  given  by 

V » F(J’)  - F(J")  (3A) 

where  J"  is  the  quantum  number  corresponding  to  the  lower  energy  state 
and  J'  is  the  quantum  number  corresponding  to  the  higher  energy  state. 

Using  Equation  32  one  obtains 

V - B(J+l)(J+2)  - BJ(J+1) 

- B<J^+3J+2)  - B(J^+J) 

- 2B(J+1)  (35) 

Since  the  lower  levels  corresponding  to  J = 0,  1,  2,...  are  generally 
populated,  absorption  of  radiation  occurs  at  frequencies 

V » 2nB,  n « 1,2,3,...  (36) 


This  equation  indicates  that  the  rotational  spectrum  observed  for  a linear 
molecule  represented  by  the  rigid  rotator  approximation  is  a set  of  lines 
spaced  by  the  factor  2B.  In  reality  the  molecules  are  not  rigid  but 
flexible  and  can  be  expected  to  stretch  under  the  influence  of  centrifugal 
forces  arising  from  their  rotation.  Thus,  the  length  of  the  molecule  and 
its  moment  of  inertia  I,  would  increase  as  the  rotational  energy  increases 
causing  B to  decrease.  The  stretching  of  the  molecule  from  a distance  r^ 
to  some  new  distance  r '-om  the  point  of  rota,  on,  gives  rise  to  a potential 
energy  of  the  form. 


V - % kx"' 


« % k(r-r^>' 


(37) 


16 


I 


! 


I 


li 


AFWAL-TR-Sl-1280 


The  rotational  energy  of  the  system  is  then  given  by  the  Hamiltonian  for 
the  system 

Ej^  = T + V 

“ % + % kCr-r^)^  (38) 

2 

Using  Equations  25  and  28,  the  first  term  becomes  P /2yr  and  Equation  38 
can  be  expressed  as 

\ (39) 

It  is  now  necessary  to  express  the  rotational  energy  E^  explicitly  in 

terms  of  the  distorted  quantity  r^v  According  to  basic  mechanics,  a 

rotating  system  of  particles  with  reduced  mass  u and  angular  velocity  w, 

2 

experiences  a centrifugal  force  urw  which  is  balanced  by  a restoring 
force  f = kr,  thus 

k(r-r^)  - 

r(k-iiw^)  - kr 

o 


Y— ”--V 

\k-yr 


(41) 


Substituting  for  r in  the  first  term  of  Equation  39,  one  obtains 

* o^nnr-ry 

o 

(1  _ _2vwi  + v?w  j ^ (k/2)(r-r  )2 
k k^  ® 


2ur^ 
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2 2 

From  (yw/k)  = (r-r^/r)  , the  last  term  in  the  first  parenthesis  can  be 
neglected  and  one  has 


2iir 


2„r 

P 


2iir. 


pV 


kr. 


J-+ (k/2)(r-r^)' 


+(k/2)(r-r^)2 

kii 

+(k/2)(r-r^)2 


ky‘r 


(42) 


From  Equation  40, 


(k/2)(r-r^)' 


, 2 4 2.2 
(u  r w ) 


2ku^r^ 


(urv^)^ 

2k 


(43) 


2 

Since  the  term  k(r-rQ)  /2  is  included  to  account  for  the  molecular 
stretching  it  is  very  small  in  comparison  to  the  kinetic  energy  term  and 
r can  be  approximated  by  r^  in  Equation  43.  Thus,  Equation  42  becomes 


r~2~ 

ky  r. 


2ku 


2pr. 


2k-.i 


(44) 
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To  convert  the  classical  result  to  a quantum  mechanical  result,  one  pro- 
ceeds as  before  using  the  rule  that  the  angular  momentum  be  quantized 
according  to 

p = /j(J+l)  h/2TT 

Substituting  for  P in  Equation  44  results  in  the  expression 


h^J(J+l) 


In  terms  of  wave  number  units 


4.  'h  T 
32Tt  kr^  u 
o 


(45) 


h 


j(j+i)  - 


32?Vk( 


Using  Equation  33  and  the  relationship 

w “ (l/2irc)  </k/m 

Equation  46  becomes, 


F(J)  - BJ(J+1)  - J^(J+1)^ 

F(J)  - BJ(J+1)  - DJ^(J+1)^ 


(46) 


(47) 


3 2 

where  0 = 4B  /w  is  known  as  the  rotational  constant  and  is  obtainable 
from  spectral  measurements.  If  one  substitutes  for  h and  c in  Equation  32, 
B has  the  value 


. 27.9_8i.  X 
I 
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Therefore,  it  is  seen  from  Equation  47  that  the  correction  term  is  very 
small  and  for  low  values  of  J the  energy  of  rotation  remains  virtually 
unchanged..  From  Equation  34  the  frequencies  of  the  allowed  transitions 
in  wave  numbers  are 

V - F(J+1)  - F(J) 

- B[(J+l)(J+2)  - J(J+1))  - 

D[(J+l)^(J+2)^  - 

- 2B(J+1)  - 4D(J+1)^  (49) 

This  shows  that  the  distance  between  successive  rotational  lines  are  no 
longer  constant  and  decrease  slightly  with  increasing  values  of  J. 

b.  Asymmetric  Top  Molecules 

Molecules  which  do  not  possess  rotational  axis  higher  than 
two-fold,  are  asymmetric  tops  and  have  the  property  of  possessing  three 
different  principal  moments  of  inertia.  The  total  angular  momentum  P 
conforms  to  the  basic  rule  that  it  remains  constant  in  magnitude  and 
direction  during  the  rotational  motion.  However,  there  is  no  longer  a 
direction  in  the  molecule  along  which  the  total  angular  momentum  has  a 
constant  component  and  there  is  no  axis  fixed  to  the  molecule  that  carries 
out  a simple  rotation  about  P.  Under  these  conditions,  according  to  the 
laws  of  classical  mechanics  (or  quantum  mechanics),  the  energy  levels  of 
an  asymmetric  top  molecule  such  as  H^O  cannot  be  represented  by  explicit 
formulation.  To  arrive  at  a set  of  qualitative  energy  levels,  one  may 
treat  the  asymmetric  top  molecule  in  terms  of  two  limiting  cases  of  the 
much  simpler  symmetric  top  molecule.  In  this  molecule,  two  of  the  prin- 
cipal moments  of  inertia  are  equal  and  denoted  by  Ig  and  the  third  by  I^.; 

For  an  asymmetric  top  molecule,  I,  ^ Ir,  where  I.  > !„  > I_.  By 

Hot  A D L 

using  the  limiting  case  in  which  Ig  is  allowed  to  decrease  gradually  from 
Ig  to  I„,  and  then  to  increase  gradually  to  I^,  Herzberg  has  shown  that 
one  can  obtain  at  best  a rough  approximation  to  the  energy  levels  of  an 
asymmetric  top  molecule.  One  of  the  more  refined  treatments  for  the 
energy  formulation  has  been  worked  out  by  Wang  (Reference  7),  For  the 
energy  term,  Wang  gives 

F(J^.)  = %(B+C)  J(J+1)  + |A  - %(B+C)]W^  (50) 
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where  the  subscript  t is  used  to  indicate  that  there  are  2J+1  energy 
levels  for  each  value  of  J.  is  a quantity  interrelated  in  a very 
complex  way  to  the  rotational  constants  A,B,C  and  the  rotational  quantum 
number  J.  Algebraic  equations  for  the  solution  for  have  been  formulated 
and  these  can  be  found  on  pp,  46-47  of  Reference  6. 

3.  ROTATION-VIBRATION  SPECTRA 

Pure  vibration  or  pure  rotation  spectra  for  polyatomic  molecules  is 
seldom  found  in  nature,  particularly  in  the  Earth's  lower  atmosphere. 

The  spectra  which  occurs  is  the  result  of  the  simultaneous  excitation 
of  rotation  and  vibration  commonly  referred  to  as  rotation-vibration 
spectra,.  For  complicated  molecules  such  as  H^O,  no  closed  formulae 
exists  as  was  explained  in  the  previous  section.  However,  mathematical 
relationships  have  been  formulated  to  fit  observed  experimental  results. 
This  is  particularly  true  of  the  H^O  molecule  whose  spectra  has  been 
measured  extensively  (References  3-5). 

As  a first  approximation,  the  vibration-rotation  energy  term  of  a 
molecule  is  given  by  the  sum  of  the  individual  vibrational  and  rotational 
terms , 

1 = G(v)  + F(J)  (51) 

where  G(v)  and  F(J)  are  given  by  Equations  22  and  50,  respectively.  The 
effect  of  the  interaction  can  be  accounted  for  by  introducing  the  effects 
of  vibration  into  the  rotational  constants  A and  D.  The  rotational  cof’ 
stants  are  affected  because  the  moment  of  inertia  I changes  during  a 
vibration  such  that  the  average  value  of  I/J  in  Equation  33  is  not  exactly 
the  same  as  its  value  in  the  equilibrium  position.  For  a linear  molecule, 
the  rotational  constant  B takes  the  form 
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where  [v]  denotes  all  the  vibrational  quantum  numbers  . and  e 

denotes  the  quantity  at  is  equilibrium  value,  a.  is  a constant  that 

accounts  for  the  deviation  of  r from  r , and  the  Coriolis  interaction  of 

e 

rotation  and  vibration.  The  influence  of  the  vibration  term  due  to  the 
action  of  the  centrifugal  force  was  accounted  for  in  the  rotational 
constant  D introduced  earlier  in  Equation  47.  The  influence  of  D will 
be  different  for  the  different  vibrational  levels  and  this  is  accounted 
for  by  substituting  D[v]  for  D.  Therefore,  the  rotational -vibrational 
energy  for  a linear  molecule  is  given  by 

T - w^(VjL+  ~)  + S S ~)(Vj^+  % + 

1 2 i k ^2  *^2 

J(J+1)  - D[v]J^(J+l)^  (53) 

For  an  asymmetric  top  molecule,  the  rotational  constants  A,  B,  and 
C take  the  form 

A|''i  ■ *6  - E “i* 

i 

B[vl  - - E 

i 

Civi  - - E “i'" 

1 

The  complete  rotational -vibrational  term  is  then  given  by 

T » G(v)  + %(B[v]  + C{v])J(J+l)  + 

(A[v]  - %(B[vl  + C[v])]W^  (55) 

where 


(v^  + %) 

(Vj^  + %) 

(V^  + %)  (54) 


A 


^ Stt^cI  ® 
A 


B — , C * 
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G(v)  is  given  by  Equation  22,  and  e denotes  the  equilibrium  value  of  the 
quantity. 

The  observation  of  infrared  spectra  must  always  be  connected  with 

an  alteration  of  the,  dipole  moment  of  the  molecule  undergoing  a transition. 

This  can  be  readily  observed  by  taking  a look  at  the  matrix  component 

M.  . of  the  dipole  moment  given  by 
1 > J 

j ■=  J fi'I'jdv  (57) 

where  and  are  the  wave  functions  for  the  indicated  ri.jlecular 
states,  is  the  dipole  moment  and  the  integration  is  over  the  entire 
configuration  space.  If  is  a constant,  then  Equation  11  becomes 

Ri  J - *5  (58) 

rom  the  orthogonality  conditions 

[ 4<\'l'jdv  - 0 (i  j)  (59) 

Equation  58  becomes 

M ^ 4»*^\|»jdv  ■ 0 (60) 

which  says  that  if  the  dipole  moment  M is  constant,  then  the  matrix  com- 
ponent of  the  dipole  is  identically  zero.  If  this  is  true,  then  the  line 

intensity  S.  . (see  Section  III)  will  be  zero  and  no  spectra  can  be 
1 1 J 

observed. 

A vibration  defines  the  coarse  structure  of  the  bands  and  a rotation 
defines  the  fine  structure  encountered.  This  can  be  seen  from  Equations 
53  and  55.  The  term  G(v)  gives  the  band  center  and  each  vibracional  level 
is  seen  to  be  accompanied  by  a series  of  rotational  levels  conrespondi ng 
to  the  2J+1  values  of  J in  the  equation.  This  is  illustrated  at  the  top 
of  Figure  2.. 
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In  general,  the  interaction  of  vibration  and  rotation  is  not  too 
large  and  the  selection  rules  for  the  rotation-vibration  spectrum  are 
the  same  as  for  the  individual  vibration  and  rotation  spectrum.  Because 
of  this  interaction,  the  rotational  selection  rules  are  governed  by  the 
direction  of  change  of  the  dipole  moment  during  the  vibration.  A rigorous 
treatment  of  the  selection  rules  involves  usage  of  the  symmetry  properties 
and  the  rotational  eigenfunction  for  each  class  of  molecules.  Such  a 
treatment  has  been  performed  for  both  the  CO2  and  molecules  by 
Dennison  (Reference  5)  also  see  Herzberg  (Reference  6 pp.  50-59  and 
pp.  468-487).  A rotational  level  is  called  positive  (+)  or  negative  (-) 
depending  on  whether  the  rotational  eigenfunction  remains  unaltered 
(symmetric)  or  changes  its  sign  (antisymmetric)  by  a reflection  of  all 
sets  of  particles  at  the  origin.  By  convention,  the  symbol  -t-*-  is  used  to 
indicate  allowed  combinations  and  -'/*■  to  indicate  the  contrary.  Transitions 
between  levels  have  to  obey  the  selection  rules 

a:  AJ*0,  ±1,  J-OV^J^O 

b:  +•*"*■  - , + •*/♦+,  - V*  “ 
c:  I • 0,  ±1 

Transitions  for  AJ  =0,  +1,  and  -1  correspond  to  Q,  R,  and  P branches, 
respectively.  The  top  of  Figure  2 shows  a simplified  illustration  of  a 
vibrational-rotational  transition  band  along  with  tie  appearance  of  the 
P,  Q,  and  R branches.;  The  Q branch  lies  at  the  band  center  with  the  P 
branch  above  the  band  center,  and  the  R branch  below  the  band  center. 

a.  CO^  Molecule 

Due  to  its  linear  symmetry,  the  CO^  molecule  has  four  funda- 
mental modes  of  vibration,  two  of  which  are  degenerate.  This  is  shown 
in  Figure  3. 

The  first  mode  corresponding  to  is  classified  as  a symmetric 
stretching  since  the  two  oxygen  atoms  move  away  symmetrically  from  the 
carbon  atom..  Such  a motion  does  not  result  in  a change  of  the  dipole 
moriient  and  is  not  infrared  active..  The  second  motion  resulting  in  is 
classified  as  a bending  mode  and  is  infrared  active..  However,  it  is  also 

1 
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Figure  3.  Fundamental  Modes  of  Vibration  of  the  CO^  Molecule 


I 

i 


a degenerate  mode  since  the  molecule  can  vibrate  with  the  same  frequency 
in  a plane  perpendicular  to  the  plane  represented  by  the  solid  arrows. 
The  third  motion  resulting  in  is  classified  as  a symmetric  stretch 
which  results  in  a change  of  the  dipole  moment  and  is  infrared  active. 

By  convention,  the  asymmetric  and  symmetric  stretching  vibrations  are 
described  as  parallel  vibrations,  where  the  change  in  dipole  moment  is 
parallel  to  the  axis  of  symmet'^y  of  the  molecule.  The  bending  vibration 
is  described  as  a perpendicular  vibration  since  tne  change  of  dipole 
moment  during  the  vibration  is  perpendicular  to  the  major  axis  of 
symmetry. 

The  infrared  spectrum  of  the  C0«  molecule  consists  principally  of 
two  intense  bands,  one  centering  at  667.38  cm  and  the  other  at  2349.15 
cm'\  In  addition  to  these  strong  bands  there  exists  a large  number  of 
weaker  bands  which  are  the  result  of  overtones  of  the  fundamentals  and 
combinations  of  overtones  and  fundamentals.  It  is  the  overtone 

2v2  (v  » 1388.18  cin‘^) 
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which  combines  with  to  give  the  combination  band  . This  band 

is  centered  at  960.97  cm”^  (10.41  pm)  and  contains  the  CO,  laser  lines 

12  16  ^ 
of  interest  for  the  C 0^  ° isotope. 

b.  H^O  Molecule 

The  H^O  molecule  has  three  fundamental  modes  of  vibration  which 
are  illustrated  in  Figure  4.  The  vibrational  energy  levels  belong  to  one 
of  two  symmetry  classes,  depending  on  whether  the  upper  energy  state  v^ 
is  odd  or  even.  Because  of  the  closeness  of  and  vibration 

resonances  do  occur,  but  they  are  of  different  symmetry  class  and  cannot 
interact  directly;  consequently,  interaction  is  very  weak.  Infrared 
spectra  has  only  been  observed  for  and  v^,  thus  the  strong  water  vapor 
bands  have  centers  at  1595  cm  *,  3755.8  cm*l , and  several  of  the  overtones 
and  combinations. 


H 


/ \ 
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Figure  4.  Fundamental  Modes  of  Vibration  of  the  H^O  Molecule 
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As  was  previously  discussed,  the  asymmetric  top  molecule  is  charac- 
terized by  the  three  unequal  moments  of  inertia  I.,  I.  where  I.  < In  < I-- 
The  three  principal  axis  about  which  the  moments  1^,  and  act  are 
denoted  a,b,c,  respectively.  The  rotational  eigenfunction  is  a function 
of  the  orientation  of  this  system  of  axis  with  respect  to  a fixed  coordinate 
system.  A condition  imposed  is  that  i/t  must  remain  unchanged  or  only  change 

K a K r 

sign  for  such  a rotation.  The  rotations  are  labeled  where 

the  subscript  2 denotes  that  the  molecule  belongs  to  the  point  group  having 
a two-fold  axis  of  symmetry  and  the  superscript  indicates  the  principal 
axis  for  the  rotation  operation.  Accordingly,  molecules  belonging  to  this 
point  group  have  three  types  of  band  contours  called  A,  B,  or  C bands 
depending  on  whether  the  change  in  dipole  moment  is  along  the  direction 
of  least,  intermediate,  or  largest  principal  moments  of  inertia.  The  H^O 
molecule  has  its  permanent  dipole  moment  along  the  intermediate  axis  of 
the  principal  moment  of  inertia  giving  rise  to  spectra  of  type  B bands. 

This  band  is  characterized  by  the  lack  of  a strong  central  branch,  with 
the  gap  filled  by  weak  lines  and  surrounded  on  both  sides  by  strong 
groups  of  lines 

For  the  frequencies  of  interest  to  CO^  laser  propagation,  the 
vibration  bands  do  not  contain  the  frequencies  of  CO^  laser  emission. 

The  extinction  due  to  H^O  is  directly  attributable  (except  for  the  R 20 
line  at  10.25  pm)  to  the  absorption  of  the  water  vapor  continuum  around 
10  pm.  This  continuum  is  attributed  to  be  the  result  of  the  accumulated 
wings  of  the  water  vapor  lines  of  the  6.3  pm  fundamental  band  and  from 
the  pure  rotational  lines  above  14  pm. 

4.  SPECTRA  PERTURBATIONS 

There  are  two  important  factors  which  have  an  effect  on  the  observed 
spectra.  These  are  termed  Coriolis  interaction  and  Fermi  resonance.  The 
degree  of  interaction  is  a function  of  the  frequency  of  the  interacting 
levels,  with  the  interaction  being  more  pronounced  for  energy  levels  lying 
closely  together. 
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The  Coriolis  interaction  comes  from  the  Coriolis  force  experienced 
by  a rotating  particle.  When  a Coriolis  force  component  happens  to  line-up 
with  the  direction  of  vibration  of  the  nuclei  in  the  molecule,  Coriolis 
interaction  will  result  with  some  coupling  of  energy  between  the  modes. 

The  closer  in  frequency  the  vibrating  modes,  the  higher  the  interaction 
between  the  energy  levels.- 

The  Coriolis  effect  is  accounted  for  in  the  energy  equation  by  the 
constant  a found  in  the  rotational  constant  B of  Equation  52,  where  a is 
the  sum  of  the  harmonic  and  anharmcrn'c  components,  and  of  the  Coriolis 
interaction.  The  interaction  is  then  accounted  for  by  replacing  a in 
Equation  52  with 

a - * (61) 

It  is  seen  from  the  energy  equation  that  the  Coriolis  interaction  in- 
creases quadratically  with  J,  such  that  for  high  values  of  J,  the  pertur- 
bation of  the  vibration  levels  become  more  significant  where  the  frequency 
of  the  levels  involved  are  spaced  closely  together.  An  example  of  the 
significance  of  this  interaction  is  given  by  Meredith,  et  al ..  (Reference  8) 
in  which  they  give  results  for  line  strengths  measured  in  the  10580  cm"^ 

and  10570  cm  ^ region.  They  found  some  line  strengths  which  are  25  times 
stronger  than  predicted  by  line  model  calculations. 

Fermi  resonance  may  occur  when  two  vibrational  levels  belonging  to 
different  vibrations  (or  combination  of  vibrations)  have  nearly  the  same 
energy.  The  occurrence  of  such  a resonance  interacts  to  cause  a pertur- 
bation of  the  energy  levels.  For  example,  such  a perturbation  arises  in 
the  CO^  molecule  for  the  state  = 1337  cm'^  and  2v2  = 1334  cm  \ The 
interaction  causes  one  of  the  levels  to  be  shifted  up  and  the  other  down 
so  that  the  separation  of  the  levels  is  larger  than  predicted,  and 
Equation  53  no  longer  allows  for  an  accurate  prediction  of  the  energy  in 
this  case.. 
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For  the  rotating-vibrating  molecule,  the  Fermi  resonance  will,  in 
general,  produce  a change  in  the  effective  rotational  constants  B of  the 
two  interacting  levels.  Correction  to  the  energy  term  can  be  accomplished 
by  introducing  a corrertion  to  the  rotational  constant.  Adel  and  Dennison 
(Reference  9)  have  formulated  such  a correction.  In  terms  of  the  unper- 
turbed rotational  constants  B°  and  B°  for  the  levels  v.  and  v.,  the 

Vi  V.  1 j’ 

effective  rotational  constants'^are 

« a^B°  + b^B® 

"'j  ""i 

B - (62) 

i j i 

2 2 

where  a and  b are  fractional  constants  obeying  the  rule  a + b =1.  A 

different  technique  in  use  to  account  for  the  Fermi  interaction  is  given 

in  Reference  10.  This  is  based  on  the  addition  of  the  extra  term 
3 3 

(J+1)  to  the  energy  equation,  where  is  a constant  determined 
from  the  rotational  constants.; 
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SECTION  III 
MOLECULAR  EXTINCTION 


1 INTRODUCTION 

The  amount  of  radiation  propagated  between  any  two  points  in  the 
atmosphere  is  determined  by  the  two  princiole  types  of  constituents, 
namely  gaseous  molecules  and  aerosols.  For  a homogeneous  optical  path 
of  length  R having  an  attenuation  coefficient  for  monochromatic 
radiation  of  wavelength  x,  the  atmospheric  transmittance  is  given  by 

T » exp(-  (Tr)  (63) 

X 

where  the  attenuation  coefficient  a is  the  sum  of  the  absorption  (k) 

A 

and  scattering  («)  coefficients  of  the  molecules  and  aerosols  encountered 
in  the  propagation  path.  Thus,  the  extinction  coefficient  is  the  sum  of 
the  four  components, 

“m 

where  the  subscript  m refers  to  molecular  and  a refers  to  aerosols.  Each 
of  these  coefficients  is  wavelength  dependent.  In  addition,  the  scattering 
coefficient  is  highly  dependent  on  the  size  of  the  scatterer  in  relation 
to  the  wavelength.  For  the  CO2  infrared  wavelengths  of  interest,  the 
molecular  contribution  to  the  scattering  coefficient  is  of  no  consequence 
and  may  be  neglected. 

This  section  deals  with  the  aspects  of  molecular  extinction  only. 

The  aspects  of  aerosol  extinction  are  treated  in  the  following  section.. 

2.  MOLECULAR  ABSORPTION 

Molecular  absorption  of  the  CO^  laser  wavelengths  is  primarily 
caused  by  the  CO^  and  H^O  molecules.  At  higher  elevations,  the  0^ 
molecule  with  a band  centered  around  9.7  pm  becomes  important.  CO^  laser 
propagation  is  significantly  affected  by  the  atmospheric  CO^  absorption 
band  which  contains  the  same  frequencies  produced  by  CO^  laser  emission. 
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Thus,  there  is  the  unusual  situation  that  the  frequencies  of  a properly 
adjusted  laser  are  exactly  coincident  with  the  center  of  the  absorption 
lines.  Absorption  due  to  H^O  is  caused  by  the  accumulated  wings  of  water 
vapor  lines  of  the  6.3  ym  fundamental  band  and  the  pure  rotational  lines 
in  the  far  infrared. 


The  molecular  contribution  to  absorption  of  radiation  by  the  atmosphere 
is  most  readily  described  in  terms  of  the  molecular  extinction  coefficient 
k^(v)  which  gives  the  fractional  energy  lost  per  unit  length  in  the 
direction  of  propagation.  The  extinction  coefficient  depends  not  only  on 
the  frequency  but  also  on  the  temperature  and  pressure  of  the  atmospheric 
constituents  along  the  propagation  path.  To  investigate  the  functional 
dependence  of  on  v,  T,  and  P,  it  is  convenient  to  define  the  absorp- 

tion coefficient  as  the  product  of  the  line  shape  factor  f(v'-v  ) and  the 
1 i ne  strength  S as 

n 

■ 2 (65) 

1-1 


where  the  summation  is  used  to  indicate  that  the  determination  of  the 
absorption  coefficient  k at  some  frequency  v requires  the  contributions 
from  several  absorption  lines  f^.  with  line  strengths  . The  line  strength 
provides  a measure  of  the  transition  probability  at  the  center  frequency. 

It  provides  a quantum  mechanical  description  giving  the  probability  that 
a photon  of  frequency  will  be  absorbed  at  the  center  frequency  of 
the  particular  transition 


E"^  + h 


oi 


(66) 


The  line  shape  f(v-v^)  relates  the  probability  S^.  of  the  transition  at 
line  center  to  that  of  any  other  frequency  v'  which  is  removed  from 
1 ine  center.. 
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For  a transition  from  state  i to  state  j,  the  line  strength  is 
given  by  (Reference  11) 


N. 


*ij 


N 


3hc 


(67) 


where 


N 


j 


N 


8i 


the  molecular  concentration  in  the  lower 
state 

total  molecular  concentration 
the  statistical  weight 
transition  frequency 
absolute  temperature 

matrix  component  of  the  dipole  moment  • 


The  temperature  and  pressure  dependence  of  the  line  strength  is  given 
by  (Reference  12) 


S - (^)“  exp  ^1  ' (68) 

° T k T T 

o 

where  S is  the  line  strength  at  standard  temperature  and  pressure  (STP) 
® -2  -1 

in  units  of  cm  atm  , P is  the  absorber  pressure  in  atm,  and  is  the 
reference  temperature..  The  value  of  m depends  on  the  molecular  species. 


3.  LINE  SHAPES 

In  examining  transitions  between  discrete  molecular  energy  levels, 
it  was  assumed  that  the  resulting  spectra  consisted  of  individual  mono- 
chromatic lines.  In  reality,  this  is  not  the  case  and  the  line  shape 
is  dependent  on  the  conditions  under  which  absorption  (or  emission)  by 
molecules  occur.  In  the  Earth's  atmosphere,  the  spectral  profile  is 
influenced  by  the  processes  of  natural  broadening,  collision  broadening, 
and  Doppler  effect.  Natural  broadening  sets  the  minimum  width  of  the 


33 


AFWAL-TR-81-1280 


spectral  line  and,  as  will  be  shown,  it  can  be  neglected  in  comparison 
to  the  other  two  processes. 


a.  Natural  Broadening 


Natural  linewidth  is  produced  by  the  interaction  of  the  molecule 
with  a radiation  field.  Transition  of  a molecule  from  one  quantum  energy 
state  to  another  occurs  over  a finite  time  interval.  It  follows  from 
Heisenberg's  uncertainty  principle,  that  the  vibration  process  accompanying 
this  transition  represents  an  infinite  set  of  frequencies  over  a specific 
spectral  range.  The  width  of  the  interval  depends  on  the  time  At  that  the 
molecule  remains  in  the  excited  state  E.  The  uncertainty  in  the  energy 
is  given  by 


h 

AE  2 

2TrAt 


(69) 


with  the  incremental  energy  involved  having  the  value 

AE  “ hcA\)  (70) 

Comparing  Equations  69  and  70  one  gets  the  result 

Av„  - crn‘^  (71) 

^ 2iicAt 

Using  the  accepted  value  of  0.1  second  for  the  lifetime  for  vibrational 

energy  states  in  the  mid-infrared, Equation  71  yields  a value  of  Av^  ' 

-11-1  ^ 

5x10  cm  which  is  insignificant  when  one  considers  CO^  laser  fre- 
quencies around  944  cm'^  . 

b.  Collision  Broadening 

This  type  of  interaction  is  also  referred  to  as  pressure 
broadening.  It  is  the  result  of  the  interaction  between  an  absorbing 
molecule  with  neighboring  molecules  resulting  in  deformation  of  the 
energy  levels  with  the  consequent  spectral  line  broadening.  The  solution 
to  this  problem  requires  the  treatment  of  very  large  numbers  of  mutually 
interacting  systems  all  of  which  must  be  treated  simultaneously.  A 
rigorous  solution  to  this  problem  does  not  exist.  Suppose  that  while  a 
molecule  is  in  the  process  of  absorbing  or  emitting  radiation,  a collision 
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with  another  molecule  takes  place  perturbing  the  radiation  process.  After 
the  collision  which  is  assumed  to  take  place  in  a time  that  is  much  less 
than  the  time  required  to  emit  the  radiation,  the  radiation  process  pro- 
ceeds as  before  the  collision.  The  result  is  that  instead  of  observing 
a single  simple  harmonic  wave  train  that  is  infinite  in  duration,  one 
observes  a wave  train  that  is  modulated  in  frequency  and  amplitude.  Thus, 
if  one  calculates  the  intensity  as  a function  of  frequency  for  many  wave 
trains  which  have  been  interrupted  by  collisions,  it  is  found  that  the 
frequency  of  the  emitted  radiation  is  spread  out  over  a small  region 
centered  at  the  frequency  of  the  original  uninterrupted  wave  train.  A 
solution  to  the  collision  broadening  problem  is  given  by  the  Lorentz 
line  profile  (Reference  12) 


f(v-v  ) “ — 

o'  TI 


a. 


(v-v  + a\ 


(72) 


where  is  the  frequency  at  the  center  of  the  spectral  line,  and  is 
the  Lorentz  half-width  of  the  spectral  line  which  is  given  by  (Reference  12) 


a 


Lo 


i-^ii 


o 


(73) 


P and  T are  the  effective  pressure  and  temperature  of  the  specimen,  and 
n is  a constant  which  is  different  for  each  of  the  molecular  species 
involved,  For  the  CO^  molecule  n has  a value  of  0.58  and  for  the  H^O 
molecule  n has  a value  of  0.62.  The  subscript  o is  used  to  designate 
the  value  of  the  quantity  at  standard  temperature  and  pressure.. 


c.  Doppler  Broadening 

Doppler  broadening  is  produced  by  the  relative  motion  of  the 
molecules  while  undergoing  an  absorption  (or  emission)  process.  Consider 
all  the  molecules  in  a gas  which  are  absorbing  radiation  during  a given 
transition.-  Because  of  thermal  motions  of  the  molecules,  some  will  be 
moving  toward  the  observer's  reference  frame  and  some  will  be  moving 
away.  The  radiation  received  from  molecules  moving  av«y  from  the  observer 
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will  appear  to  be  of  lower  frequency.  Conversely,  the  molecules  moving 
toward  the  observer  will  appear  to  emit  higher  frequencies.  The  result 
is  that  the  apparent  width  of  the  emission  is  broadened.  The  shape  of 
a Doppler  broadened  line  is  given  by  (Reference  11) 


f(v-v^)  « (£n2/Tr)^(l/ajj)exp 


-(v-Vo)^An2 


(74) 


where  is  the  Doppler  half-width  given  by 

- 3.58  X 10“^v^(T/M)^  (75) 

T is  the  temperature  at  which  the  process  takes  place,  and  M is  the 
molecular  weight  of  the  gas. 


At  sea  level,  pressure  broadening  attains  a maximum  and  decreases 
with  altitude  in  proportion  to  the  reduction  in  pressure  encountered. 
Near  sea  level  the  Doppler  line  width  is  negligible  in  comparison  to  the 
half-width  arising  from  pressure  broadening.  At  some  altitude  the  width 
of  a pressure  broadened  line  will  equal  the  Doppler  width.  Above  this 
point,  the  Doppler  width  becomes  increasingly  larger  and  is  the  dominant 
mechanism  for  line  broadening.  For  altitudes  above  15  kilometers,  both 
Doppler  and  pressure  broadening  have  to  be  accounted  for  simultaneously.; 
This  leads  to  the  Voight  line  profile  (Reference  12)  given  by 


f(v) 


“1 


[(V-V^) 


2opX  2 2 

J + ® T 

ln2  ^ 


(76) 


4.  WING  THEORY 


A comparison  of  the  Doppler  and  Lorentz  line  profiles  with  similar 
strengths  and  line  widths  is  shown  in  Figure  5.;  An  important  observation 
is  that  the  Gaussian  shape  of  the  Doppler  profile  allows  only  a small 
absorption  in  the  wings  and  thus  the  wings  of  a molecular  absorption  line 
is  mainly  the  result  of  Lorentz  broadening..  Since  the  Voight  profile  is 
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a convolution  of  the  two  curves  shown,  this  means  that  the  resulting  line 
profile  should  be  dispersed.  The  degree  of  exactness  in  regards  as  to 
how  well  the  far  wings  are  described  by  this  dispersion  profile  has  not 
been  resolved,  and  line  wing  theory  awaits  more  experimental  data  before 
an  exact  model  can  find  general  acceptance..  As  a result,  several  shape 
models  have  been  proposed..  The  most  prominent  one  is  the  Van  Vleck- 
Weisskopf  model  given  by  Reference  13 


f(v-v^)’ 


-)I 


“L 


“l 


-]  (77) 


where  v is  the  resonant  frequency  of  the  molecule  and  is  the  half-width 

of  the  line  at  half-maximum  intensity.:  In  the  infrared,  absorption  is 

significant  only  in  the  region  <<  Thus,  in  this  region 

(\j+v  ) > {v-\>  ) and  \)/v  ' 1.  Under  these  conditions  Equation  77  reduces 

0 0 0 

to  Equation  72. 


a 


Figure  5,  Doppler  and  Lorentz  Spectral  Line  Shapes  for 
Similar  Intensities  and  Line  Widths 
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The  water  vapor  continuum  has  been  and  still  is  a source  of  contro- 
versy. The  most  accepted  theory  is  that  the  continuum  absorption  is  due 
to  the  influence  of  the  remote  wings  of  the  6.3  pm  absorption  band  and 
the  rotational  lines  above  14  pm  of  the  water  molecule..  Recent  workers 
(References  14,  15)  have  hypothesized  a second  absorption  mechanism  in 
the  continuum  and  have  advanced  the  dimer  (H20)2  theory.  Measurements 
by  Arefev  (Reference  16)  of  the  attenuation  of  CO^  laser  radiation  by 
water  vapor  at  different  temperatures  and  concentrations  point  in  favor 
of  the  dimer  theory.  The  dimer  theory  is  still  not  accepted  by  the 
community  and  continuum  models  in  ucs  ore  based  solely  on  contributions 
from  the  wings  of  the  water  vapor  lines.  Several  workers  have  investi- 
gated the  contributions  from  the  far  wings  in  the  continuum  (References 
17-20),. 

A study  by  Roberts,  Selby,  and  Bieberman  (Reference  20)  of 
attenuation  in  the  8-14  pm  region,  has  led  them  to  postulate  the 
empirically  derived  line  shape  factor 


(78) 


where  a is  a constant  obeying  the  rule  l<a<2...  Roberts  et  al . have 
shown  that  Equation  78  models  the  contribution  from  the  line  wings  much 
more  precisely  than  the  Lorentz  model. 


5,  TRANSMISSION  MODELS 

Several  atmospheric  transmission  models  have  evolved  during  the 
yea'^s.  Most  of  these  are  based  on  the  band  model  concept  in  which  com- 
puted spectral  lines  are  averaged  into  spectral  bands.  A spectral  band 
is  assumed  to  contain  an  array  of  spectral  lines  of  a certain  shape, 
intensity,  half-width,  and  spacing  according  to  some  simple  statistical 
model.  In  real  spectra,  this  is  not  the  case  and  this  is  the  major  source 
of  discrepancy  between  real  and  computed  spectra. 
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For  the  computation  of  laser  atmospheric  extinction,  line  models 
are  used..  The  extinction  coefficient  is  given  by  Equation  65  where 
f(v-v^)  is  the  line  shape  factor.;  For  a single  line  where  the  Lorentz 
line  profile  is  applicable,  the  extinction  is  given  by 

S ctj 

k(v)  - 1— (79) 

TT  “L 

In  general  a large  number  of  absorption  lines  belonging  to  different 
molecules  contribute  to  the  attenuation  at  any  specific  laser  frequency. 
In  this  case  the  total  optical  depth  (O.D.)  requires  evaluation  and  is 
given  by 

O.D.  “ k(v)W 


i vv  ^i.j  °i.j  “j 

Ti  ^ 4—*  2 2 

j-i  1-1  “ i.j’ 


(80) 


where  m^  is  the  amount  of  the  j-th  absorbing  gas,  W is  the  total  absorber 

amount  expressed  in  reciprocal  units,  and  S.  . is  the  line  strength  given 

1 » J 

in  Equations  67  and  68.  The  double  summation  accounts  for  the  fact  that 
there  are  m separate  absorbing  species  each  having  n lines  with  significant 
absorption  at  frequency  v.  This  is  the  form  of  equation  on  which  the  high 
resolution  transmission  model  HITRAN  (Reference  19)  is  based.;  Clough 
(Reference  21)  et  al  . have  developed  an  accelerated  version  of  HITRAN 
called  FASCJDE,  which  cuts  the  computer  time  required  by  a factor  of 
approximately  100.  They  have  achieved  this  by  developing  an  algoritnm 
whicn  accelerates  the  convolution  of  the  line  shape  function  with  the 
spectral  line  data  used.  The  contribution  of  continuum  absorption  for 
self  and  foreign  broadened  components  are  included  in  the  model.  A 
layered  atmosphere  is  assumed.  The  Voight  line  profile  is  utilized  for 
all  layers  with  a bound  of  25  cm”^  except  at  the  higher  altitudes  where 
the  Doppler  line  profile  is  used.  The  accuracy  for  the  model  is  given  as 
nominally  ±0.5“'.  except  in  the  wings  of  Voight  lines  for  which  the  error 
may  reach  ±37;.  FASCODE  does  well  in  calculating  the  contribution  due  to 
molecular  extinction,  but  as  of  this  writing,  this  model  does  not  account 
for  aerosol  attenuation  in  the  middle  infrared.; 
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SECTION  IV 
AEROSOL  EXTINCTION 


1 . INTRODUCTION 

The  attenuation  of  radiation  by  scattering  is  caused  by  discrete 
variations  in  the  refractive  index  due  to  the  presence  of  particles  in 
the  medium.  The  main  scatterers  are  water  droplets,  dust,  smoke,  and 
haze  particles.  If  the  incident  radiation  has  a wavelength  much  smaller 
than  the  size  of  the  scattering  particles,  the  interaction  of  the  radiation 
with  the  scatterers  results  in  Rayleigh  scattering.  If  the  dimensions  of 
the  scattering  particles  are  comparable  to  wavelength  of  the  incident 
t^adiation,  the  interaction  cf  the  radiation  with  the  scatterers  results 
in  Mie  scattering.  At  the  infrared  wavelengths  of  interest,  Rayleigh 
scattering  is  not  applicable  in  describing  the  process  and  one  has  to  use 
Mie  scattering  theory.;  The  exact  solution  to  the  Mie  scattering  problem 
is  extremely  difficult  with  the  main  difficulty  being  that  one  cannot 
express  the  angular  scattering  by  a simple  analytical  expression,  the 
phase  function  can  only  be  described  as  a series  of  Legendre  polynomials 
and  this  makes  the  numerical  solution  very  laborious.  Basically,  Mie 
theory  relates  the  scattering  and  absorption  of  electromagnetic  energy 
by  spherical  particles  having  complex  index  of  refraction  m = n^  * jn^.; 

From  Mie  theory,  Oeirmendjian  (Reference  22)  has  shown  that  the  extinction 
efficiency  factorsQ(A)  are  given  by 

00 

Qs  - (2n+l)<|a  + |b„|2)  (81) 

n“l  “ ” 

00 

+ V (82) 

n-1 

“ ^ext  * (8.'^) 
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where  the  subscripts  a and  s refer  to  absorption  and  scattering  respectively, 
and  is  the  extinction  (or  attenuation)  efficiency  of  the  attenuating 
particles.  The  quantity  x is  shown  as  the  size  parameter 

X = 2TTr/X  (84) 

which  is  a function  of  both  the  wavelength  and  radius  of  the  particle.; 

The  quantities  a and  b are  coe  ficients  in  the  expression  for  the 
scattered  and  attenuated  fields.  The  quantities  a^  and  b can  be  expressed 
in  terms  of  spherical  Bessel  functions  and  Hankel  functions  of  the  second 
kind  with  arguments  x and  complex  refractive  index  m, 

- <j'j^(B))j>'^(x) 

where  e -■  mx  and  <t)  = i|;  + ja.. 

The  quantities  i|)  and  n are  Ricatti-Bessel  functions  and  are  defined  in 
terms  of  half  integral  order  Bessel  and  Neumann  function  with  argument  x, 

♦«(*>  - <-7-)'* 

n(x)  - 

The  extinction  efficiency  factors  as  given  by  the  exact  Mie  theory  are 
difficult  to  compute  and  even  with  modern  computers,  the  effort  is  time 
consuming.;  For  computational  purposes,  the  approximations  worked  out  by 
Van  de  Hulst  (Reference  23)  and  modified  by  Deirmendjian  (Reference  24) 
are  much  more  convenient  to  use.  Van  de  Hulst  obtained  the  following 
formulas  for  the  extinction  and  absorption  efficiency  factors. 


a_ 


n 


n 
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'^ext  ■ 2 - 4e''’*^“®(5^)8ln(o-6) 


-4e-^‘“B(S^)2co.(p-26) 


+4(£^)2  cos  28 


(85) 


*abs 


--4xn, 

1 + e 1 

2xno 


-4xn 


1-1 


8x^2 


(86) 


where 


tan  6 


4Trr 

A 

"2 

n^-1 


(n^-l) 


2Trr 


and  and  are  the  real  and  imaginary  part  of  the  refractive  index, 
respectively.  The  above  expressions  were  derived  with  the  assumptions 
(n-1)  <<  1,  <<  1,  and  Zitr/x  >>1 . Because  of  the  assumptions  in  his 

theory,  Van  de  Hulst  approximations  overestimate  the  extinction  cross 
section  for  .mall  p and  underestimate  it  in  varying  degrees  as  p approaches 
and  surpasses  the  position  for  the  first  maximum  of  at  p = 4,08. 
Deirmendjian  (Reference  23)  worked  out  correction  factors  which  when 
applied  Lc  the  Van  de  Hulst  equations  give  a good  approximation  to  the 
exact  Mie  equations,  Deirmendjian  (Reference  24)  has  shown  that  the  use 
of  his  correction  factors  D cn  the  Van  de  Hulst  equations  give  values  that 
are  within  6^^  of  the  values  obtained  from  Mie  calculations  as  applied  to 
water  particles  at  wavelength  near  10  pm.  In  view  of  the  uncertainty  in 
dtop-size  distributions  and  sphericity  of  water  particles,  the  corrected 
Van  de  Hulst  equations  are  widely  used  and  will  be  used  in  this  study. 
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However,  the  Dermendjian  correction  factors  have  been  modified  to  give 
better  approximations  for  wavelengths  near  10.6  ym.  The  Van  de  Hulst 
equations  as  modified  by  Deirmendjian  give  larger  errors  for  particles 
with  radius  smaller  than  1.0  ym  (48%  for  r = 0.5  ym).  The  modified 
Deirmendjian  correction  factors  give  answers  which  are  within  1.0%  of 
the  exact  Mie  theory.;  Deirmendjian's  correction  factors  and  the  improved 
Deirmendjian  factors  for  wavelengths  near  10.6  ym  are  discussed  in  Appendix 
A For  the  case  of  fogs  and  precipitation  at  the  10.6  ym  wavelength  of 
interest,  the  most  recent  values  of  n^  = 1.185  and  n^  = 0.0662  given  by 
Hale  and  Querry  (Reference  25)  will  be  used  as  opposed  to  Centenno's 
(Reference  26)  values  of  n^  = l.;18  and  = 0.08, 

The  power  AP  removed  from  a laser  beam  by  a particle  with  cross- 
sectional  area  vr^  is  given  by 

AP  - PQ(X.r)Trr^  (87) 

For  a unit  volume  containing  N particles  of  the  same  size  and  composition, 
the  total  power  removed  from  the  beam  becomes 

AP  - PNO(X.r)Ttr^  (88) 

In  dealing  with  atmospheric  phenomena,  the  particles  encountered  are  of 
varying  sizes  and  composition.  Assuming  that  the  particles  are  of  known 
concentration  n(r)  and  of  the  same  composition.  Equation  88  becomes 

N 

Cf'ext  ' ^ n(r^)Q(X.rj)r^^  (89) 

i»l 

where  the  attenuation  ratio  aP/P  is  known  as  the  extinction  coefficient.. 
Denoting  the  largest  and  smallest  particle  radii  found  in  the  distribution 
by  R^  and  R.|  respectively,  the  integral  representation  of  Equation  89  is 
given  by 

^ext  * ^ J r2^n(r^)Q(X,rj)dr^  (90) 
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Q(x,r.)  is  a dimensionless  quantity,  therefore  the  units  of  are 
determined  by  the  units  of  the  radius  r and  the  volume  concentration  n(r). 

_3 

For  r in  ym  and  n(r.)  in  cm  , will  have  the  more  standard  units  of 

*1  1 0Xl  2 

km”  when  Equations  89  or  90  are  multiplied  by  10"“^. 

2.  AEROSOL  PROPERTIES 

Extinction  of  laser  radiation  by  aerosols  is  the  result  of  both 
absorption  and  scattering.  To  derive  accurate  results  from  models,  one 
must  know  the  density  and  distribution  of  the  size  droplets  in  the  pro- 
pagation path.  Aerosol  composition  varies  mostly  with  time  of  the  year, 
climate,  altitude,  and  geographical  location..  For  ocean  front  and  other 
marine  type  atmospheres,  salt  spray  predominates.  Over  continental  areas, 
aerosol  particles  consisting  of  a mixture  of  soluble  and  insoluble  materials 
occur.  Formations  of  mixed  particles  by  coagulation  of  aerosols  produced 
by  industrial  processes  do  occur.  The  properties  of  these  mixed  particles 
are  subject  to  change  with  changing  atmospheric  conditions.  A particle 
can  change  size  or  form  due  to  a condensation  or  evaporation  process; 
thus,  it  is  difficult  to  formulate  a model  that  holds  for  all  conditions. 
Many  analytical  expressions  have  been  formulated  to  represent  actual 
particle  size  distributions.  These  analytical  expressions  have  given  way 
to  empirically  derived  formulation  based  on  knowledge  of  the  composition 
and  properties  of  the  aerosols. 


The  effect  that  an  aerosol  particle  has  on  che  propagation  of 
radiation  in  the  atmosphere  depends  in  part  on  its  complex  index  of 
refraction 


m - n,  - jn. 


(91) 


This  is  particularly  true  for  CO2  laser  radiation  propagating  in  high 
water  vapor  concentration  where  absorption  by  water  particles  is  determined 
by  the  imaginary  part  of  the  refractive  index  n^. 
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The  physical  composition  of  continental  aerosols  varies  between  the 
extremes  of  dry  insoluble  dust  particle  (arid  soil)  and  a clear  salt 
soluble  particle  (industrial  process).  Junge  (Reference  27)  has  shown 
that  these  particles  grow  with  humidity.  Measurements  (Reference  27) 
indicate  that  these  mixed  pafticles  behave  more  like  droplets  of  salt 
solutions  with  expected  refractive  indicies  between  that  of  NaCl  and 
that  of  H^O.  The  refractive  index  of  dry  NaCl  is  in  the  vicinity  of 
1.5  while  the  real  part  of  the  refractive  index  of  H^O  is  approximately 
1,2  for  wavelengths  near  10  ym. 

Because  raindrops  are  not  particles  suspended  in  the  atmosphere, 
they  are  not  considered  aerosols..  However,  they  exhibit  atmospheric 
optical  properties  to  radiation  similar  to  those  exhibited  by  fog  par- 
ticles. For  this  study,  rain  drops  (precipitation)  will  be  treated  as 
an  aerosol . 


3.  SIZE  DISTRIBUTION 

a.  Continental  Aerosols 

One  of  the  earliest  formulations  proposed  to  represent  continental 
aerosols  is  attributed  to  Junge  (Reference  27),,  This  takes  the  form  of  a 
power  law  and  may  oe  represented  by 


n(R)  - R* ' (92) 

2,3 

where  R is  the  particle  radius,  C is  a constant,  and  v is  a parameter  with 
value  between  2.5  and  4,0,  This  distribution  characterizes  large  particles 
satisfactorily  but  it  is  not  very  good  in  characterizing  smaller  particles 
and  fails  in  the  representation  of  maritime  aerosol  distributions.  The 
distribution  also  has  a singularity  as  R->0.  An  alternate  distribution 
which  avoids  this  singularity  is  the  regularized  (RPL)  power  law  (Reference 
28) 

N(r)  - (r/a)'')  (93) 
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where  N(r)  = number  of  particles  with  radii  greater  than  r,  is  the 

total  number  of  particles  per  unit  volume,  and  a and  v are  adjustable 

parameters . 

Another  expression 

on  the  gamma  function. 

takes  the  form 

that  adequately  describes  wet  aerosols  is  based 
The  distribution  given  by  Zuev  (Reference  29) 

fCa) 

■1 

1 y+1  a'^  .-ya/r 

r(„+i) 

where 

r (y+l) 

m 

gamma  function 

y 

- 

parameter 

r 

m 

most  probable  radius 

a 

- 

particle  radius  . 

Ward,  et  al . (Reference  30)  examined  a number  of  models  anc  included 
that  the  RPL  model  provided  the  best  fit  to  the  examined  dat  Including 
their  own. 

b.  Fog  Particles 

Experimental  results  have  shown  that  the  power  law  does  not 
describe  adequately  fog  particle  distributions.  The  gamma  function  may 
be  used  to  describe  fog  particle  distributions;  however,  a simpler  dis- 
tribution formulated  by  Shifrin  (Reference  31)  has  found  wider  acceptance. 
This  is  given  by 

n(R)  - ' (95) 

where  n(R)  is  the  volume  concentration  at  the  radius  (R),  and  A,  B,  a, 
and  Y are  positive  constants.  The  parameters  a and  y can  be  adjusted  to 
fit  an  observed  fog  particle  size  distribution.  Having  chosen  a and  y. 
the  remaining  constants  can  be  uniquely  determined  by  the  total  number 
of  particles  N per  unit  volume  and  by  the  mode  radius  r^  where  the  con- 
centration is  a maximum.,  Chu  and  Hogg  (Reference  32)  modified  slightly 
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Shifrin's  distribution  by  introducing  the  term  R 
to  obtain 


r/r  in  Equation  95 
m 


1 

n(r/r„)  » A(r/r„)“ 


The  constant  A Cfin  he  determined  from  the  integral 


aj 

N » n 


(R)dR  . 


where  Equation  85  has  been  substituted  for  n(R).  Letting  y = BR^,  then 

R - (y/B)^/“ 


dR  - 


Substitution  for  R and  dR,  and  regrouping  terms.  Equation  98  becomes  the 
standard  gamma  integral 


,a+l/Y 


r(a+l/Y) 


.a+l/Y 


Thus  the  value  of  A is  found  to  be 


,a+l/Y 


r(a+l/Y5 


(100) 


The  value  of  B is  determined  by  taking  the  derivative  of  n(R)  and  setting 


it  to  zero 


dn(R) 


ArO-I  e"®RY(a-YBR’)  - 0 
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Thus,  a = yBR^  and  using  the  fact  that  R=1  for  r=r^  gives 

B - (101) 

Y 

This  particle  size  distribution  has  been  found  to  give  a good  represen- 
tation for  fog  particles  and  other  aerosols  (References  31,  32).  By  proper 
Fitting  of  the  parameters  a and  y to  the  corresponding  fog  particle  size 
distribution,  the  observed  fog  can  be  properly  modelled.  Figure  6 
illustrates  several  distributions  for  selected  values  a and  y. 

c.  Precipitation 

Precipitation  models  developed  have  been  based  on  empirical 
results  (References  32-34).  Laws  and  Parsons  (Reference  34)  found 
experimentally  that  the  rain  drop  size  distribution  is  a function  of  the 
total  precipitation  rate.  Their  results  were  not  shown  in  the  form  of  a 
distribution  function  n(a)  but  as  a related  distribution  m(a,p),  being 
given  by  the  fraction  of  the  total  volume  of  water  striking  the  ground 
due  to  drops  of  radius  a.  Goldstein  (Reference  35)  gives  the  relationship 
between  the  two  distributions  as 

n(a)  (4n/3)v(a)a^ 

mCa.p)  . ^ ^ (102) 

4Tt/3  I n(a)v(a)a-^da 

o 

where  v(a)  is  the  terminal  velocity  of  the  drops  with  radius  a in  m/sec, 
Goldstein  determined  that  the  precipitation  rate  is  given  by 

«D 

p - 15.1  n(a)v(a)a^da  (103) 

Substituting  this  expression  in  Equation  102  one  obtains 

15.1  n(a)v(a)a^ 

ni(a,p)  « (104) 

P 

where 

v(a)  = 9.6(1  - e’^1®)  (105) 
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Figure  6 


REATiVE  DROP  RADIUS  (/^m) 

Size  Distribution  for  Fog  Model 


49 


AFWAL-TR-81-1280 


and  a is  in  cm.  Using  Equation  104  to  substitute  for  n(a)  in  Equation  90, 
one  obtains 

<D 

f Q£(^,a)  m(a,p)da 


f QE(X.a) 

- piT  I 

j 15.1 


v(a)a 


Defining 


g(a,X)  - 


wQE<X,a) 

15.1  v(a)a 


(106) 


the  extinction  is  given  by 


A 

■ i 


g(a,X)  ia(a,p)da  (km  ) >•  (107) 


Observations  by  Marshall  and  Palmer  (Reference  36)  led  to  the 
Marshall-Palmer  distribution 


(108) 


Nn  - N e"^0  * 
u o 

where  u • number  of  drops  in  a xinlt  volume  with 


diameters  in  the  interval  D and  D + dD  > 
empirically  derived  constant  with 
value  of  0.08  cm“^, 
constant  related  to  rainfall  rate,  K 
in  units  of  mm/hr  and  having  a value 
of  A.lR'^’^^cm”^. 
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A comparison  of  the  Marshal  1 -Palmer  and  the  Lav/s  and  Parson  distributions 
are  plotted  in  Figure  7 versus  experimental  observations  for  three  different 
rain  rates.  The  plots  indicate  that  the  Palmer-Marshall  distribution 
does  not  fit  observations  for  drops  of  small  diameters  and  that  the  deviation 
gets  larger  with  increasing  rain  rates.  However,  it  should  be  noted  that 
as  the  rain  rates  increase  the  drop  size  also  increases  and  the  indicated 
difference  in  the  distribution  will  be  minimized.  While  the  Laws  and 
Parsons  formula  gives  a better  fit  to  experimental  results  for  smaller 
sizes,  it  requires  knowledge  of  the  additional  quantity  v^.  - The  terminal 
velocity  v^  can  be  approximated  by  (Reference  37) 

" 1420  cm/ sec  (109) 

for  drops  with  diameters  between  1 and  4 mm.  Chen  (Reference  38)  formulated 
Laws  and  Parsons  relationship  as 

n.  - (110) 

‘ v.v. 


where  v^.  is  the  volume  of  the  particle.  For  the  range  of  validity  of 
Equation  109,  Equation  110  can  be  approximated  by 


•2.5 


n. 


8410 


(111) 


It  is  important  to  remember  that  these  distributions  are  empirical 
equations  based  on  accumulation  of  highly  variable  samples.  By  examination 
of  rainfall  data  obtained  during  a two-year  period,  Geotis  (Reference  39) 
found  that  at  diameters  greater  than  0.6  mm  (the  smallest  size  he  could 
measure),  the  drop  concentration  was  exponentially  distributed.  He  reported 
that  the  exponential  in  Palmer-Marshall  distribution  varied  with  rainfall 
in  an  unpredictable  manner.  He  found  that  A = 3.7  gave  a good  fit  to 
the  observed  dropsize  distribution  curve,  where  is  defined  as  that  drop 
diameter  which  divides  the  drop  distribution  in  such  a way  that  half  the 
water  content  is  contained  in  drops  greater  than  D^.  The  drawback  is  that 
cannot  be  measured  directly. 
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Figure  7.  Comparison  of  the  Marshall-Palmer  Distribution 
(Solid  Lines)  and  Observed  Results  of  Laws  and 
Parsons  (Broken  Lines)  for  Rain  Rates  of  1 , 5, 
and  10  mm/hr.  Adapted  from  Battan,  Reference 
37,  p.  41 

d.  Snow 

Snow  is  a form  of  precipitation  which  is  difficult  to  characterize 
because  of  the  irregular  shape  of  the  sn.’^w  flakes.  A form  of  the  Marshal  1- 
Palmer  distribution  given  by  Equation  108  has  been  used  to  characterize 
snow.  The  value  of  the  parameter  D refers  to  the  diameter  of  a sphere  of 
water  having  the  same  mass  as  the  original  particle.  When  using  the 
applicable  quantities  (Reference  40) 

N - 0.025R’®*9^  cm"^ 


AFWAL-TR~81-1280 

in  Equation  108,  the  size  distribution  for  snow  particles  is  given  by 

Njj  « 0.025R"0*9^  exp(2.29R“®*^^)  (115) 

where  R is  expressed  in  units  mm"\ 
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SECTION  V 

SYSTEM  DESIGN  AND  ANALYSIS 


1,.,  INTRODUCTION 

Measurements  performed  in  an  uncontrolled  environment  such  as  the 
atmosphere  require  special  techniques  and  instrumentation.  Where  practical, 
possible  sources  of  error  should  be  eliminated  by  suitable  design  of  the 
test  apparatus.  Where  this  is  not  possible  (or  practical),  the  design 
should  minimize  the  error  involved.  Because  the  propagation  path  used 
crosses  a 200  feet  taxiway  and  lies  alongside  a runway,  atmospheric 
turbulence  creates  problems  which  require  special  attention.  This  section 
covers  the  effects  caused  by  atmospheric  turbulence,  system  design,  and 
analysi s . 

2.  ATMOSPHERIC  TURBULENCE 

Atmospheric  turbulence  is  identified  with  fluctuations  in  air  tem- 
perature and  pressure  which  causes  fluctuations  in  the  transmitting  medium 
index  of  refraction.  Since  the  pressure  (and  humidity)  change  appreciably 
at  a given  altitude,  most  of  the  refractive  index  fluctuations  arise 
primarily  from  temperature  fluctuations  in  the  atmosphere,  which  are 
produced  by  t^e  turbulent  mixing  of  hot  air  with  cool  air..  This  mixing 
occurs  in  zones  or  packets  called  eddies  which  range  in  size  from  the 

inner  scale  of  turbulence  1 , to  the  outer  scale  of  turbulence  L . 

0 0 

The  outer  scale  of  turbulence  L^  is  a quantity  associated  with  the 

introduction  of  energy  in  the  turbulence.  It  is  defined  as  that  minimum 

size  of  large  scales  at  which  turbulent  energy  is  no  longer  introduced 

in  the  field.  The  inner  scale  of  turbulence  is  the  size  of  the  smaller 

scale  that  results  in  the  conversion  of  energy  in  the  turbulence  into 

heat.  The  outer  scale  L varies  in  size,  from  tenths  of  a meter  to  a 

0 

few  meters  near  ground  level,  while  the  inner  scale  1^  is  the  order  of 
one  millimeter  near  ground  level  (Reference  41). 
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and  are  also  referred  to  as  atmospheric  inhomogeneities . It 
is  these  inhomogeneities  with  their  varying  refractive  indexes  that  causes 
refraction  and  diffraction  of  the  optical  wavefront. 

The  index  of  refraction  for  the  turbulent  layer  can  be  written  as 
the  mean  of  the  index  of  refraction  n in  air,  and  a fluctuation  An  about 
the  meaf, 

Hj,  - 1 + An  (116) 

where  the  refractive  index  of  air  has  been  rounded  to  n = 1 . The  fluc- 
tuations are  correlated  only  over  short  distances  and  hence  produce  a 
general  degradation  in  optical  signals  propagating  through  the  atmosphere. 
In  more  specific  terms,  the  fluctuations  form  turbulent  eddies  which 
cause  refraction  and  diffraction  of  the  optical  wavefront.;  To  see  this 
more  clearly,  consider  a plane  optical  wavefront  propagating  i-'  air  between 
a receiver  and  a transmitter  as  shown  in  Figure  8.  Let  the  wavefront  be 
incident  at  an  angle  with  respect  to  a plane  turbulent  layer  of  air  of 
thickness  L whose  refracMve  index  is  given  by  Equation  116. 


Figure  8,  Projection  of  a Diffracted  Laser  Beam  Due  to 
Turbulent  Eddies 
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After  the  wavefront  enters  the  turbulent  eddy  L,  diffraction  causes  the 
beam  to  be  deflected  at  an  angle  A0.  away  from  the  intended  trajectory. 
The  new  angle  of  propagation  after  entering  the  turbulent  eddy  is 


e 


i 


+ Ae. 

O 1 


The  coherence  length  L is  defined  as  the  eddy  size  over  which  the 

index  of  refraction  is  approximately  constant.  By  the  definition  of  the 

coherence  length,  the  turbulent  eddy  of  size  in  Figure  8 acts  as  a 

dielectric  wedge  with  respect  to  the  optical  wavefront  rather  than  as  a 

plate.  Therefore,  the  direction  of  the  optical  wavefront  after  emerging 

from  the  turbulent  eddy  is  identical  to  the  direction  within  the  eddy 

which  differs  in  direction  from  that  of  incident  radiation.  The  effects 

repeat  themselves  for  subsequent  eddies  along  the  path  so  that  the  overall 

effect  is  a summation  over  the  different  coherence  intervals  encountered. 

Since  each  angular  deviation  is  uncorrelated  with  the  angular  deviation 

in  the  previous  coherence  interval,  the  total  mean  squared  angular 
2 

deviation  <A0  > is  a function  of  the  sum  of  the  individual  mean  squared 
angular  deviations. 


When  a laser  beam  passes  through  a turbulent  atmosphere,  its  inter- 
action  with  the  eddies  produces  random  variations  in  amplitude  and  phase 
of  the  signal.  This  randomness  leads  to  a variety  of  effects  which  tend 
to  degrade  the  performance  of  a laser  system.  Although  somewhat  inter- 
related, the  most  dominant  effects  are  beam  steering,  spot  (image)  dancing, 
beam  spreading,  image  blurring,  and  scintillation. 

a.  Spot  Dancing 

When  the  laser  beam  passes  through  pockets  of  air  (eddies)  in 
which  the  refractive  index  is  different  from  the  norm,  the  entire  beam 
is  deviated  from  its  original  direction.  Depending  on  the  frequency  and 
number  of  these  eddies  along  the  path,  the  result  is  a movement  of  the 
laser  beam  in  the  receiver  plane,  which  in  turn  causes  the  image  to  be 
focused  at  different  points  in  the  focal  plane.  Chiba  (Reference  h?) 
gives  a relation  for  calculating  spot  dar:ing  (also  known  as  beam  wander). 
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According  to  Chiba,  the  standard  deviation  of  the  linear  displacement 
of  the  laser  beam  in  the  receiver  plane  is  given  by 

(117) 

where  is  the  beam  size  at  the  transmitter  output,  X is  the  laser  wave- 
length, L is  the  distance  between  transmitter  and  receiver,  and  is  the 
index  structure  constant  which  can  be  obtained  from  a measurement  of  the 
temperature  structure  constant  C^,  using  the  relationship  (Reference  43) 

(118) 

This  quantity  does  not  lend  itself  to  accurate  measurements.  Estimates 

derived  by  Tatarski  from  experimental  data,  compared  with  measurements 

from  Genoud  and  Goodwin  give  an  order  of  magnitude  agreement  (R,  erence 

44).  Davies  calculated  that  for  turbulence  near  the  ground,  can  be 

approximated  by  the  value  corresponding  to  intermediate  turbulence  and 

-8  -1/3 

gives  the  approximate  value  of  4 X 10"  m * . For  strong  turbulence 

-7  -1/3 

= 5 X 10  m ; thus,  it  can  be  seen  that  there  exists  a large 
range  of  values  for  and  computations  involving  this  quantity  are  at 
best  very  rough  estimates. 


b.  Beam  Spreading 

Interaction  of  the  laser  with  the  inhomogeneities  in  the  atmo- 
sphere results  in  small  angle  scattering.  The  net  effect  is  that  localized 
inhomogeneities  produce  changes  in  the  divergence  of  the  beam  with  a 
resulting  spread  of  the  signal  energy  over  a larger  area.  The  long-term 
spread  for  a Gaussian  beam  may  be  approximately  by  the  relation  (Reference 
45) 

®8  " 


where  k if,  the  wave  number  and  is  the  lateral  coherence  length  given  by 


P 

o 


1.74  S lO'^X 
C 

n 


(120) 
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c.  Scintillation 

Known  also  as  power  fluctuation,  scintillation  is  the  effect  wherein 
destructive  interference  within  the  beam  produced  local  fluctuations  in  the 
amplitude,  so  that  across  the  beam  there  appears  dark  and  bright  areas 
compared  to  the  average  level  of  brightness.  This  leads  to  fluctuations 
in  the  received  power  which  are  strongly  dependent  on  the  receiving  aper- 
ture dimensions.  Tatarski  has  treated  the  subject  of  power  fluctuations 
quite  extensively  in  his  book  (Reference  46). 

Tatarski  postulates  that  the  amplitude  variations  obey  a log  normal 
distribution  and  derives  the  following  expression  for  the  fluctuation  of 
the  logarithmic  level  of  the  amplitude 

- log  (A/A^)2 

R 

- C (r)Z^^^^dZ  (121) 

y ti 
o 

which  holds  when  L >>  XR  >>  £«.  A is  the  initial  amplitude,  k is  the 
0 0 0 

wavenumber,  i and  I are  the  inner  and  outer  scales  of  turbulence  and 
0 0 

R is  the  path  length.  Equation  121  holds  for  a propagation  path  of  varying 
altitudes,  and  it  allows  for  changes  in  the  turbulent  structure  along  the 
path.  For  paths  which  are  nearly  parallel  to  the  earth's  surface  along 
regions  wnere  the  turbulence  structure  does  not  change  appreciably,  Equation 
121  reduces  to 


^ - .31  <122) 

a 

The  fluctuations  in  the  logarithmic  level  of  intensity  I is  related  to 
that  of  amplitude  A by 

» 4 (123) 
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Tatarski  defines  the  maximiim  observable  scintillation  frequency  in 
terms  of  the  correlation  interval  For  the  case  of  a receiver  with 
an  aperture  d <<  R,  the  maximum  scintillation  frequency  given  by 


max 


(i24) 


where  is  the  wind  speed  normal  to  the  line  of  sight  between  the  trans- 
mitter and  receiver.  For  d >>  the  higher  spatial  frequencies  will 
be  discriminated  against,  and  the  frequency  spectrum  will  be  significant 
up  to  a frequency  approximately  expressed  by  (Reference  47) 


f - Vjd  (125) 

d.  Beam  Steering 

Beam  steering  is  the  deviation  of  the  beam  from  the  Intended 
line  of  sight.  The  lateral  displacement  of  the  beam  has  an  rms  deviation 


y - Rog/2  (126) 

along  the  axis  of  the  receiver  plane,  where  may  be  approximated  by 

Equation  119  (Reference  48) 

0 - Ti/kpQ 


Thus, 


y - Rii/2kp^ 


- RX/4p„  (127) 

e.  Image  Blurring  ” 

Image  blurring  is  defined  as  the  destruction  of  the  phase 

coherence  across  the  beam,  causing  the  phase  to  change  rapidly  with 

position  resulting  in  a blurring  of  the  image.  The  angular  blur  of  a 

2 

point  source  is  Guassian  with  zero  mean  and  fluctuation  o0  given  by 
Equation  119.; 
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3.  SYSTEM  DESIGN 

In  the  design  of  a laser  transmissometer,  there  are  three  possible 
choices  for  the  selection  of  the  receiver  aperture  size  D . The  aperture 
can  be  less,  equal,  or  greater  than  the  beam  diameter  d.  The  first  case 
(D^  < d)  requires  precise  knowledge  of  the  position  of  the  receiving 
aperture  in  relation  to  the  beam  geometric  configuration.  Because  of 
the  position  accuracies  required  to  track  a wandering  beam,  calibration 
in  terms  of  total  power  received  becomes  very  difficult  and  can  lead  to 
very  large  errors.  The  second  case  (D^  = d)  suffers  from  the  same  con- 
siderations and  measurement  errors  can  be  worse  than  the  first  case  due 
to  the  turbulence  effects  discussed  earlier.  By  selecting  > d,  most 
of  the  turbulence  effects  can  be  eliminated.  Blurring  and  beam  dancing 
effects  will  be  present  but  they  can  be  minimized  by  proper  receiver 
design. 

a.  Gaussian  Optics 

The  Gaussian  distribution  which  describes  the  fundamental  mode 
of  a laser  is  given  by  (Reference  49) 

I(r)  - (128) 

2 

iTa 

where  is  the  total  output  power,  I(r)  is  the  intensity  at  the  distance 

r from  the  beam  center  and  a is  the  radius  where  the  intensity  falls  to 
2 

1/e  or  13.5%  of  the  central  intensity.  The  diffraction  produced  by 

passing  a Gaussian  beam  through  a circular  aperture  which  truncates  the 

distribution  at  the  1/e  points  has  been  treated  by  Dickson  (Reference  50). 

In  the  absence  of  aberrations  a lens  irradiated  by  a Gaussian  beam  will 

form  an  image  that  is  dependent  on  the  ratio  of  lens  aperture  to  beam 

diameter  2a.  As  the  aperture  is  increased,  the  diffraction  effects  are 

reduced  and  as  Dickson  shows  (Reference  50),  diffraction  effects  are 

entirely  eliminated  when  the  aperture  diameter  is  4a.  Because  of  the 

high  cost  of  large  diffraction  lim’.  ?d  optics  (and  system  growth),  it 

is  the  accepted  criteria  that  a well  designed  system  may  be  truncated 
2 

at  the  1/e  intensity  points.  This  is  achieved  with  an  aperture  which 
is  at  least  1.4  times  the  beam  width. 


60 


AFWAL-TR-81-1280 


The  far  field  beam  divergence  of  a laser  beam  is  given  by 

2X  (129) 


0 


Tia 


2 . 


where  e is  the  full  angle  beam  divergence  measured  at  the  1/e  intensity 
points  and  a is  the  beam  radius  before  propagation.  Due  to  the  long 
wavelength  of  a CO^  laser,  the  beam  divergence  is  rather  large,  and  the 
beam  diameter  will  increase  accordingly  to  the  relation  d = Le,  By  in- 
creasing a before  transmission,  a corresponding  reduction  in  0 is  achieved. 
This  is  achieved  by  the  use  of  a beam  expander,  which  gives  a collimation 
range  of 

z z 

(130) 


2 2 
inn  a 


This  relation  is  easily  derived  from  the  geometry  considerations  shown 
on  Figure  9.  As  the  beam  propagates,  it  grows  in  radius  according  to 
the  relationship  (Reference  51) 

w - w tl  )^]^ 

” nw 

o 

where  w is  the  beam  radius  at  its  waist. 

0 

b.  Laser  Transmissometer 

The  general  layout  of  the  laser  transmitter  is  illustrated  on 

Figure  10  with  the  beam  expander  configuration  illustrated  on  Figure  11. 

The  distance  from  the  laser  to  the  first  element  of  the  telescope  is 

170  cm.  At  this  first  surface,  the  beam  radius  has  grown  from  2.5  to 

3.39  mm,  and  the  minimum  aperture  to  provide  truncation  at  less  than  the 
2 

1/e  points  is  9.50  mm.  Using  Equation  130,  a trade-off  analysis  between 

range  and  magnification  yields  the  selected  magnification  of  20X.  This 

establishes  the  minimum  size  of  telescope  output  aperture,  at  190  mm  to 

2 

keep  truncation  at  less  than  the  1/e  points.  Ignoring  diffraction  and 
atmospheric  effects,  use  of  Equation  130  gives  a beam  width  of  19  cm  at 
a range  of  1874  meters  and  a beam  width  of  13.5  cm  at  the  half  range 
point.;  At  the  selected  range  of  1000  meters,  the  beam  width  is  approxi- 
mately 14.7  cm.. 
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Figure  10.  Optical  Layout  of  Laser  Transmitter 
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The  laser  transmitter  (Figure  10)  consists  of  a multi-lineCO^  laser 
which  is  mounted  coaxially  to  a HeNe  laser.  A is  an  attenuator-polarizer 
combination  used  to  present  a safe  level  of  power  to  the  receiver  detector 
during  low  extinction  conditions.  The  first  component  attenuates  the 
signal  and  the  second  is  used  to  restrict  transmission  to  the  perpendicular 
component  of  the  E vector.  A two-degree  wedge  is  used  to  axially  combine 
the  CO^  and  HeNe  laser  outputs.  The  signals  are  modulated  at  a IkHz  rate 
using  a WWVB  NBS  reference  time  signal  which  is  also  used  to  phase-lock 
the  receiver  to  the  transmitter.  B^  is  a 97%  reflective  beam  splitter 
which  provides  a small  portion  of  the  signal  to  the  reference  detector 
Pm.  A 20X  beam  expander  is  used  to  provide  an  8-inch  beam  at  the  receiver. 
The  energy  at  the  receiver  is  collected  by  a 24-inch  f/5  parabolic  mirror 
which  is  tilted  slightly  off-axis.  The  signal  is  folded  and  passed  through 
a couple  of  beam  splitters  which  are  used  to  limit  the  energy  reaching  the 
detector  (Figure  11).  In  addition,  the  first  beam  splitter  is  used  to 
provide  a reference  signal  for  calibration  of  the  receiver  and  to  measure 
the  received  HeNe  whenever  desired.  Lens  L is  used  to  relay  the  CO^ 
signal  into  an  f/10  solid  recollector.  This  arrangement  provides  a steady 
signal  to  the  detector  during  periods  of  high  turbulence.  When  properly 
aligned,  signal  variation  is  approximately  1%  in  periods  of  extreme  tur- 
bulence. Pyroelectric  detectors  are  used  at  both  receiver  and  transmitter. 
This  limits  the  dynamic  range  of  the  system  to  72dBs.  The  signal  pro- 
cessing electronics  are  shown  in  Figure  11.  The  signal  conditioning  and 
recording  are  almost  identical,  except  that  in  the  receiver  section,  a 
phase-lock  amplifier  is  used  to  phase-lock  the  receiver  to  the  transmitter. 

The  IkHz  modulating  signal  was  selected  because  of  convenience  for 
phase-locking  and  because  it  is  well  beyond  the  maximum  scintillation 
frequency  that  could  be  experienced.  Using  the  most  pessimistic  case  of 
Equation  124  and  a wind  velocity  of  30  m/hr,  the  maximum  scintillation 
frequency  was  found  to  be  130  Hz. 

Weather  and  particle  size  data  was  collected  at  the  receiver  end 
only.  To  ensure  a nearly  isotropic  atmospheric  path,  the  receiver  signal 
was  visually  monitored  using  a digital  voltmeter.  Stable  voltmeter 
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readings  indicated  a good  probability  of  a stable  atmospheric  path,  while 
a varying  voltmeter  reading  indicated  an  unstable  atmospheric  path. 


There  are  several  factors  which  must  be  considered  in  selecting  an 
optimum  aperture  size.  These  are: 

1.  Optical  aberrations 

2.  Diffraction  effects 

3.  Atmospheric  effects 


The  transmitter  telescope  uses  an  off-axis  parabola  with  a spherically 
corrected  first  element.  This  combination  yields  a virtually  aberration 
free  image  of  the  magnified  beam.  The  diffraction  effects  are  eliminated 
by  choosing  an  aperture  at  least  twice  the  beam  diameter.  For  a beam 
diameter  of  14.7  cm,  the  diffractive  component  of  the  aperture  is  29.4  cm. 


As  previously  mentioned,  use  of  the  index  structure  constant  requires 
a very  good  knowledge  of  the  types  of  turbulence  encountered.  Because  the 
experiment  is  conducted  in  an  area  where  high  turbulence  is  a large  proba- 
bilit)',  the  value  of  = 5 X 10'^  will  be  used  in  this  study.  This 
gives  a coherence  length  of 


1.74  X 10‘2  X (10.6)2  X 10"^2 

[ . 5.44  cm 

10-^  X 25  X 10“!^ 


Spot  dancing  is  calculated  with  the  help  of  Equation  117.  Using 
7 1/3 

= 5 X 10'  m"  ^ in  Equation  120,  the  coherence  length  for  an  optical 
path  length  is  16.3  cm.  Spot  dancing  is  found  to  be 


(1.90 


2.85  cm 
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Using  = 5.44  cm  in  Equation  119,  the  atmospheric  spread  is  found  to  be 

« _ LX 


8 zp. 


9.74  cm 


Beam  steering  has  two  components  which  must  be  accounted  for.  The  first 
may  be  termed  a short-term  component  and  may  be  mathematically  approximated 
by  Equation  130.  This  is  found  to  be 


y - 


" 4.87  CO 

The  long-term  component  is  caused  by  changes  in  the  refractive  index  of 
the  atmosphere  due  to  temperature  variations  and  atmospheric  particulate 
content.  The  day  to  day  variations  in  the  refractive  index,  change  the 
path  traversed  by  the  beam  and  results  in  the  beam  falling  at  a different 
position  in  the  receiving  aperture.  This  effect  could  be  minimized  by 
using  a true  off-axis  parabola,  but  because  of  the  cost  involved,  this 
is  not  practical.  The  only  practical  solution  is  to  relay  all  the  energy 
captured  by  the  aperture  into  the  detector.  This  may  be  accomplished  by 
realigning  the  laser  transmitter  prior  to  taking  measurements  so  that  the 
beam  is  at  the  optical  center  of  the  aperture.  This  beam  deviation  may 
be  treated  as  an  alignment  error  of  linear  displacement  d = Re,  where  6 
is  the  refracted  angle  from  the  intended  optical  path.  Assuming  a maximum 
value  of  0 = 100  microradians,  the  increase  in  aperture  diameter  is  10  cm. 

Collecting  all  the  terms  involved  which  contribute  to  enlargement 
of  the  receiver  aperture,  the  minimum  diameter  is  found  to  be 

D - d + [p^j^  + - 41.5  cm  • (132) 


i] 


II 


f 
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The  chosen  aperture  has  a 61  cm  diameter  which  satisfies  all  of  the  above 
criteria  for  a range  of  up  to  2.2  kilometers  if  the  long-term  drift  d is 
compensated  before  initiating  measurements. 

4.  MEASUREMENT  ERRORS 

The  quantity  that  is  measured  directly  is  the  transmittance  T.  The 
extinction  is  computed  from  the  transmittance  value  T and  from  Equation  63, 

'Tj-'l/Llnl'  (133) 

where  the  path  length  L has  the  value  of  1 kilometer  and  a has  the  units 

- 1 ^ 
of  km  . To  obtain  in  dB/km,  one  has  to  multiply  the  obtained  quantity 

by  4.34.  This  is  shown  in  Appendix  B. 

The  error  involved  in  the  calculations  is  less  than  2%  and  this  is 
due  to  the  uncertainty  in  the  readout  of  the  calibration  device.  To 
minimize  this  error,  the  calibrating  radiometer  was  calibrated  each  time 
prior  to  calibration  of  the  transmissometer,  and  the  same  instrument  was 
used  to  calibrate  both  receiver  and  transmitter. 

For  moderate  to  high  extinction  conditions,  the  2%  error  in  trans- 
mittance becomes  negligible  when  computing  extinction.  This  is  easily 
seen  by  expressing  the  natural  log  of  T in  a power  series  and  then  taking 
the  total  differential  of  o . 

A 

- dll  - %(T  - 1)  + R.O.T]  dT 

V 

- -dT  + 0(T  - 1)”  ' (134) 

Retaining  the  linear  term  only  and  dividing  by  o , one  obtains  the  approx- 

A 

i mat ion 

d^M  2 1 dT 

"rf"  I 

InT 

If  one  assumes  a 2X  error  in  the  measurement  for  T = .01,  the  error  in 

extinction  is  in  the  order  of  0.4%  while  for  T = lO”^  the  error  is  in  the 

-4 

order  of  0.14%.  For  the  high  extinction  cases  where  T > 10  an  error  of 
10®  gives  =>rrors  in  extinction  in  the  order  of  1%.. 
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SECTION  VI 

DISCUSSION  OF  RESULTS 


1.  INTRODUCTION 

Results  of  extinction  measurements  are  presented  with  measurements 
supported  with  particle  size  distributions  where  applicable.  For  clear 
and  summer  hazy  weather,  the  measured  transmission  (or  extinction)  values 
are  compared  with  those  given  by  the  FASCODE  model.  The  values  obtained 
with  the  model  are  of  sufficient  accuracy  to  determine  the  true  molecular 
extinction,  but  as  pointed  out  in  Section  III,  this  model  cannot  account 
for  aerosol  extinction.  Clear  weather  results  are  presented  first, 
followed  in  sequence  by  fog,  rain,  haze,  and  snow. 

Measurements  were  made  at  discrete  times  at  single  wavelengths  and 
during  very  short  periods  over  several  lines.  With  the  piezoelectric 
mounted  mirror,  it  was  possible  to  make  transmission  measurements  over 
six  different  wavelengths  in  approximately  60  to  75  seconds.  This  short 
time  was  necessary  to  obtain  extinction  measurements  at  several  wavelengths 
under  the  same  (or  nearly  the  same)  atmospheric  conditions. 

2.  CLEAR  WEATHER 

Three  examplesof  clear  weather  measurements  are  tabulated  in  Tables 
1,  2,  and  3.  Data  in  Table  1 is  representative  of  a humid,  cold  winter 
day,  while  the  data  in  Table  2 is  representative  of  a dry  winter  day. 

Both  measured  and  calculated  values  indicate  that  for  low  temperatures, 
transmission  is  not  greatly  affected.  Table  3 shows  measured  and  calcu- 
lated transmission  values  for  a hot,  humid  summer  day.  For  the  cold 
winter  days,  the  transmission  losses  are  dominated  by  the  CO2  contribution 
to  extinction.  Fo<"  the  hot  humid  days,  extinction  is  dominated  by  the 
water  vapor  continuum.  Two  points  to  be  observed  are  that  for  the  cold 
winter  days  the  measured  transmission  is  approximately  2%  lower  than  that 
calculated  using  '^ASCODE,  and  for  the  hot,  humid,  summer  days  the  measured 
transmission  is  slightly  higher  than  that  calculated  by  the  FASCODE  model 
The  2%  difference  in  transmission  values  for  the  winter  measurements  may 
in  reality  be  less  than  that.  Calculations  by  McCoy  (Reference  52)  give 
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TABLE  1 

COMPARISON  OF  MEASURED  AND  CALCULATED 
TRANSMISSION  FOR  CLEAR  WEATHER  CONDITIONS, 
25  FEBRUARY  1981,  1300  HRS.  VISIBILITY  WAS 
APPROXIMATELY  18  KM 


LINE 

NO. 

A (vim) 

Tjjj 

(measured) 

Tc 

(computed)  AT 

P28 

10.S75 

.914 

.925 

.011 

P26 

10.653 

.909 

.922 

.013 

P22 

10. tu 

.900 

.910 

.010 

P20 

10,591 

.895 

.906 

.011 

P18 

10.571 

.897 

.906 

.009 

R18 

10.260 

.893 

.906 

.013 

R16 

10.275 

.888 

.900 

.012 

Temperature 

Dew  Point  ■ 

- 36*F 

34’F 

Relative  Humidity 

Pressure  “ 990  mb 

- 82:: 

I 

I 
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TABLE  2 

COMPARISON  OF  MEASURED  AND  CALCULATED 
TRANSMISSION  FOR  CLEAR  WEATHER  CONDITIONS, 
3 MARCH  1981,  0830  HRS.  VISIBILITY  WAS 
GREATER  THAN  25  KM 


LINE  NO. 

X (ym) 

^m 

(measured) 

Tc 

(computed) 

AT 

P28 

10.675 

.549 

.969 

.020 

P26 

10.653 

.945 

.965 

.020 

P22 

10.611 

.935 

.957 

.022 

P20 

10.591 

.931 

.953 

.022 

P18 

10.571 

.930 

.952 

.022 

P14 

10.532 

.927 

.951 

.024 

R18 

10.260 

.926 

.94i7 

.021 

R16 

10.275 

.922 

.945 

.023 

Temperature  " 28*F  Relative  Humidity  • 40% 

Dew  Point  ■ 17 ’F  Pressux'e  - 990  mb 
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TABLE  3 

COMPARISON  OF  MEASURED  AND  CALCULATED 
TRANSMISSION  FOR  CLEAR,  HUMID,  SUMMER 
WEATHER  CONDITIONS,  20  JULY  1981,  1300  HRS, 
VISIBILITY  WAS  GREATER  THAN  20  KM 


LINE  NO. 

A (urn) 

(measured) 

(computed) 

AT 

P28 

10.675 

.588 

.572 

.016 

P26 

10,653 

.584 

.570 

.014 

P22 

10.611 

.575 

.560 

.013 

P20 

10.591 

.572 

.556 

.016 

P18 

10.571 

.574 

.560 

.014 

R16 

10.275 

.576 

.562 

.014 

R18 

10.260 

.593 

.575 

.018 

Temperature  - 23.73  C Relative  Humidity  » 93.75Z 

Dew  Point  ■ 22.66  C Absolute  Humidity  - 20i.03 

gr/m3 
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ar,  absorption  of  0.06  ktn"^  for  a CO^  concentration  of  335  parts  per  million 
for  January  type  weather.  This  turns  out  to  be  a 94.2%  transmission  which 
falls  halfway  in  between  the  measured  and  computed  values  by  FASCODE. 

3.  FOG 

Over  150  fog  distributions  were  measured  during  the  course  of  this 
study.  A very  large  number  of  these  were  in  the  light  fog  category  and 
of  little  interest  to  this  study.  Size  distributions  for  three  of  these 
fogs  are  shown  in  Figures  12,  13,  and  14.  Measured  and  calculated  trans- 
mission values  corresponding  to  the  distribution  of  Figure  12  are  tabulated 
in  Table  4.  Transmission  data  for  size  distribution  of  Figure  13  is  tab- 
ulated in  Table  5,  and  that  of  Figure  14  is  tabulated  in  Table  6.  The 
fog  of  20  January  was  very  stable  and  the  measured  distribution  hardly 
changed  at  all.  The  dotted  curve  corresponds  to  a size  distribution 
measurement  taken  25  minutes  later.  The  extinction  obtained  from  Mie 
theory  yields  only  a 0.02  km"^  extinction  due  to  the  fog  particulates. 

This  is  about  a 2%  loss  in  transmission.  From  Table  5 it  is  seen  that 
the  difference  betweeen  the  measured  and  computed  values  is  approximately 
7%;  thus,  there  is  a 5%  unaccounted  difference.  A similar  observation 
can  be  made  for  the  data  obtained  on  23  January.  The  contribution  to 
extinction  by  fog  particulates  was  found  to  be  0.008  km”\  The  unaccounted 
transmission  difference  is  approximately  3%. 

Stable,  light  fog  existed  all  day  during  February  17th.  Computation 
of  the  aerosol  contribution  to  extinction  at  10.6  ym  yielded  a value  of 
0.1  km’\  This  compares  well  with  the  difference  between  the  measured 
and  computed  extinction  of  0.114  (Table  6)  at  10.6  ym. 

Heavy  fog  conditions  were  encountered  during  the  18th  and  19th  of 
February  1981.  During  those  two  days  the  extinction  exceeded  the  dynamic 
range  of  the  system  more  than  half  the  time  during  which  measurements 
were  taken.  To  help  in  the  analysis,  it  was  useful  to  apply  the  relation- 
ship between  extinction  and  the  liquid  water  content  in  a fog.  Chylek 
(Reference  53)  derived  a linear  relationship  of  the  form 

T • bW  (136) 
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TABLE  4 

MEASURED  AND  CALCULATED  TRANSMISSION  FOR 
SEVERAL  COg  LASER  LINES  UNDER  LIGHT  FOG 

CONDITIONS,  20  JANUARY  1981,  1510  HRS.  DATA 
CORRESPONDS  TO  FOG  DISTRIBUTION  OF  FIGURE  12 


LINE  NO. 

X (ym) 

Tjh 

(measured) 

Tc 

(computed) 

AT 

P28 

10.675 

.888 

.955 

.067 

P26 

10.653 

.878 

.945 

.067 

P20 

10.691 

.866 

.935 

.069 

P18 

10.671 

.864 

.932 

.068 

P14 

10.532 

.861 

.930 

.069 

RI8 

10.260 

.867 

/\^A 
• U/ V 

Temperature  • 36 ’F 

Dew  Point  ■ 29 'F 

Relative  riiraidity  ■ 407. 

Visibility  - 4 V. 
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TABLE  5 

MEASURED  AND  CALCULATED  TRANSMISSION  FOR 
SEVERAL  CO2  LASER  LINES,  23  JANUARY  1981, 

0920  HRS.  LIGHT  STABLE  FOG  WAS  PRESENT 
DATA  CORRESPONDS  TO  FOG  DISTRIBUTION  OF 
FIGURE  13 


LINE  NO. 

X (ym) 

Tin 

(measured) 

T 

■‘’C 

(computed) 

AT 

P28 

10.675 

.892 

.938 

.046 

P26 

10.653 

.889 

.934 

.045 

P22 

10.611 

.883 

.yz5 

# A 

. UHZ. 

P20 

10.591 

.880 

.920 

.040 

P18 

10.571 

OOA 
• OUU 

.920 

.040 

R16 

10.275 

.873 

.911 

.038 

RIO 

10.260 

.877 

.916 

.039 

Temperature  = 32.2®] 

Dew  Point  ■ 26. 5® F 

F Relative  Humidity  * 

Visibility  * 2 km 

757o 
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TABLE  6 

MEASURED  AND  COMPUTED  EXTINCTION  VALUES 
FOR  SEVERAL  WAVELENGTH  UNDER  LIGHT  FOG 
CONDITIONS,  17  FEBRUARY  1981,  1440  HRS 


LINE 

NO. 

X 

km“l 

^m, 

krn-l 

km“l 

kin“l 

P28 

10.675 

.291 

.172 

.119 

P26 

10.653 

.297 

.177 

.120 

P20 

10.591 

.308 

.194 

.114 

P18 

10.571 

.304 

.198 

.119 

R18 

10.260 

.302 

.194 

.108 

Temperature  ■ 53.35*F 
Dew  Point  - 51.85T 


r.eletive  Humidity  • 967, 
Pressure  - 990  mb 
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where  b is  a constant  and  W is  the  liquid  water  content.  Chylek  showed 
this  relationship  to  be  independent  of  the  size  distribution  and  also 
that  the  relationship  holds  only  for  a definite  wavelength  which  is 
determined  by  the  size  of  the  largest  droplets  present  in  a given  layer 
of  fog  or  clouds.  Chylek  results  indicate  that  the  maximum  wavelength 
for  which  Equation  136  holds  is  11  ym,  for  a distribution  in  which  the 
larger  particles  have  maximum  radii  r^  near  14  ym.  For  this  case  the 
constant  b in  Equation  136  has  a value  of  128.  Using  the  data  provided 
and  the  derived  relationship 

(137) 

^ 2pX 

the  relationship  at  10.6  ym  gives 

Tw  “ 147  W (138) 

with  the  larger  particles  having  radii  of  13  ym.  It  is  not  critical 
that  r be  less  than  r , but  that  the  larger  sizes  do  not  contribute 

iFi 

heavily  to  either  a or  W. 


Equation  136  was  derived  empirically  from  variable  samples  and  is 
not  exact.;  However,  for  the  case  of  heavy  extinction,  one  does  not  need 
to  be  exact,  and  the  approximation  it  gives  provides  one  with  additional 
information  for  correlation  with  extinction  obtained  through  use  of  Mie 
calculations . 


The  liquid  water  content  of  a concentration  of  particles  N of  size 
r,  density  p,  and  volume  V is  given  by 

W - pVN(r)  • (139) 

For  a polydispersion  of  spherical  water  Iroplets,  the  total  liquid  water 
content  is  obtained  by  integrating  Eq*  tion  139  over  the  size  range  of 
the  particle  distribution.  Equation  139  then  takes  the  form 


W " TTp  ^ ^ N(r)r^dr‘ 


(140) 
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The  number  distribution  was  converted  to  a density  distribution  by  dividing 
N(r)  by  the  volume  sample  and  then  dividing  this  value  by  the  sampling 
time  of  the  measurement.  Numerical  integration  was  then  performed  to 
obtain  W. 

Extinction  data  for  the  morning  of  18  February  1981  is  tabulated  in 
Table  7.  Between  0605  and  0655  hours,  extinction  in  the  optical  path 
exceeded  the  dynamic  range  of  the  system.  The  extinction  computed  from 
the  size  distribution  and  the  extinction  cross  section  is  tabulated  under 

3 

0^.  The  last  two  columns  give  the  fog  liquid  water  content  in  gr/m  and 

the  extinction  obtained  by  use  of  Equation  138.  Comparison  of  a and  a 

c w 

gives  an  amazingly  good  correlation  between  the  two  quantities.  Measure- 
ments taken  at  0705  and  0715  show  a good  correlation  between  o and  o 
indicating  that  the  fog  was  very  uniform  at  these  times.  Figures  15  and 
16  show  several  size  distribution  corresponding  to  the  data  tabulated  in 
Table  7.  Both  of  these  figures  illustrate  the  bimodal  tendency  of  some 
of  the  distributions.  Both  of  these  figures  show  that  a small  shift  in 
the  concentration  of  large  particles  can  cause  a large  change  in  the 
measured  extinction.  Distribution  for  0610  and  0615  hours  appears  to  be 
almost  identical  in  size  but  from  Table  7,  it  can  be  seen  that  there  is 
a 31  dB/km  difference  in  extinction  between  them.  TabK  8 gives  extinction 
data  obtained  in  medium  density  fog  and  Table  9 presents  data  collected 
in  light  density  fog.  As  can  be  observed,  there  is  not  much  change  in 
extinction  between  the  difference  wavelengths.  Whatever  differences  are 
observed  are  mostly  due  to  the  dynamic  nature  of  the  fog.  Figure  17 
illustrates  the  size  distribution  for  the  light  fog  conditions  of  the  data 
shown  in  Table  9.  Table  10  gives  the  transmission  and  extinction  values 
during  the  fog  lifting  period  for  the  18th  of  February. 

February  19th  provided  an  excellent  opportunity  for  fog  measurements. 

Fog  would  roll  in,  stabilize,  and  then  move  out  to  be  followed  by  fog  with 
entirely  different  characteristics.  Data  for  a continuous  three  and  one- 
half  hour  period  is  tabulated  in  Tables  11  through  13.  Size  distributions 
representative  of  the  data  tabulated  in  Table  11  are  plotted  in  Figures  18 
and  19.  Figure  19  shows  the  stable  nature  of  the  fog  between  0545  and 
0605  hours.  The  0550  distribution  corresponds  to  an  extinction  of  66.9  dB/km.; 
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TIME 

0605 

0610 

0615 

0670 

0625 

0630 

0635 

0640 

0645 

0650 

0655 

0700 

0705 

0710 

0715 


TABLE  7 

EXTINCTION  OF  10.6  ym  CO^  LASER  RADIATION. 
DATA  TAKEN  FROM  0605  TO  0715,  18  FEB  81  . 


r 

n 

ym 

a 

dB/fm 

dB/kffl 

V 

9 

gr/m 

°w 

dB/ksk 

3.53 

— 

140.8 

.2181 

139.1 

3.52 

— 

143.4 

.2384 

152.1 

3.64 

— 

112.6 

.1863 

118.9 

3.75 

— 

82.0 

.1257 

80.2 

3.84 

— 

61.0 

.0992 

63.3 

3.65 

— 

75.4 

.1256 

80.1 

2.92 

— 

63.2 

.1064 

67.9 

3.30 

— 

73.8 

.1239 

79.0 

3.29 

— 

93.4 

.1540 

98.2 

3.40 

_ 

102.4 

.1674 

106.8 

3.20 

— 

81.1 

.1334 

85.1 

2.61 

18.83 

22.8 

.0396 

25.3 

2.18 

10.28 

10.9 

.0712 

13.5 

2.54 

16.81 

22.8 

.0386 

24.6 

2.26 

10.45 

10.4 

.0201 

12.8 
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TABLE  8 

EXTINCTION  OF  CO^  LASER  RADIATION  IN  MEDIUM 

DENSITY  FOG  FOR  SEVERAL  LINES,  18  FEBRUARY 
1981,  0658-0659  HOURS 


X(vm) 

(dB/kn) 

^ (dB/km) 

P28 

10.675 

33.3 

P26 

10.673 

33.8 

P20 

10.591 

34.0 

P18 

10.571 

33.9 

RIB 

10.260 

35.6 

R16 

10.275 

35.6 

Temperature  - 55 ’F  Relative  Humidity  • 100% 

Dew  Point  ■ 55*F  Pressure  ■ 990  mb 
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TABLE  9 

EXTINCTION  OF  CO^  LASER  RADIATION  IN  LIGHT 

DENSITY  FOG  FOR  SEVERAL  WAVELENGTHS, 

18  FEBRUARY  1981,  0715-0716  HOURS 

LINE  NO.  X (yiQ)  ^ (dB/km) 


P28 

10.675 

8.98 

P26 

10.653 

9.28 

P20 

10.591 

10.28 

P18 

10.571 

9.41 

R18 

10.260 

10.44 

Temperature  « 54 'F 
Dew  Point  ■ 54*F 


Relative  Humidity  ■ iOOZ 
Pressure  • 990  mb 
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a 


f ' 


TABLE  10 


MEASURED 
FOG  LIFT 

TRANSMISSION  AND  EXTINCTION  DURING 

ING  PERIOD,  18  FEBRUARY  1981 

TIME 

T (10"^) 

^ (db/km) 

0717 

16.14 

7.92 

0718 

26.70 

5.74 

0719 

33.70 

4.73 

0720 

41.41 

3.83 

0721 

51.74 

2.86 

0722 

55.94 

2.58 

0723 

64.56 

1.90 

0724 

69.15 

1.60 

0725 

74.02 

1.31 

0726 

76.52 

1.16 

0727 

76.20 

1.18 

Temperature  - 53.8*F 
Dew  Point  • 54* F 


Relative  H\imidity  - 99% 
Pressure  * 990  mb 


.1 

I 


I. 


i 

i 


I 

j 

I 
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TIME 

0455 

0500 

0505 

0510 

0515 

0520 

0525 

0530 

0535 

0540 

0545 

0550 

0555 

0600 

0605 


TABLE  11 

EXTINCTION  OF  CO2  LASER  RADIATION  IN 
HEAVY  FOG.  DATA  TAKEN  19  FEB  81 


3.06 

2.59 

2.99 

2.66 

2.58 

2.15 

2.35 

3.66 

3.44 

3.81 

3.99 
3.89 

4.99 
4.59 
4.75 


a 

dB/^ 

0 

dB/&B 

dB/ln 

74.1 

.1231 

78.5 

31.5 

20.4 

.0531 

33.9 

52.7 

54.2 

.0891 

56.8 

36.9 

29.2 

.0505 

32.2 

28.3 

27.3 

.0498 

31.8 

17.2 

7.0 

.0154 

9.8 

15.1 

4.6 

.0087 

5.5 

44.3 

53.7 

.0840 

53.6 

31.7 

34.8 

.0547 

34.9 

76.7 

66.4 

66.7 

148.4 
103.1 
114.8 


.1176 

.1015 

.1025 

.2200 

.1539 

.1707 


75.0 

64.8 

65.4 

140.4 

98.2 

108.9 
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EXTINCTION  OF  CO,  LA 
FOG,  13  FEB  81  ^ 


TIME 

r 

n 

pm 

0 

dB/kn 

0610 

2.81 

34.7 

0615 

4.74 

— 

0620 

3.55 

— 

0625 

3.18 

— 

0630 

3.41 

— 

0635 

3.37 

53.9 

0640 

3.51 

52.6 

0645 

3.94 

52.9 

0650 

3.82 

51.6 

0655 

3.60 

53.2 

0700 

3.11 

38.4 

0705 

3.18 

37.9 

0710 


3.10 


37.6 


12 


RADIATION  IN  HEAVY 


0 

dB/to 

1 

1 

w 

0 

dS/km 

31.5 

.0548 

35 

64.1 

.09,  ' 

61.1 

55.8 

.0865 

55.2 

49.9 

.0795 

50.7 

64.4 

.1001 

63.9 

54.7 

.0852 

54.4 

51.9 

.0802 

51.2 

54.0 

.0828 

52.8 

49.4 

.0754 

48.1 

54.5 

.0841 

53.7 

40.6 

.0641 

40.9 

38.7 

.0613 

39.1 

36.6 

.0581 

37.1 

i 
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TABLE  13 

EXTINCTION  OF  10.6  ym  CO,  LASER  RADIATION 
IN  HEAVY  FOG,  19  FEB  81 


TCIE 

r 

m 

yn 

°in 

dB/kn 

dB/km 

V 

gr/m^ 

dB/kn 

0715 

4.25 

— 

107.4 

.1704 

108.7 

0720 

4.57 

— 

173.3 

.2749 

173.4 

0725 

4.46 

— 

167.1 

.2614 

166.8 

0730 

5.50 

— 

270.9 

.4357 

277.9 

0735 

5.66 

— 

291.0 

.4707 

300.2 

0740 

5.73 

— 

281.0 

.4566 

291.3 

0745 

5.23 

— 

276.5 

.4455 

284.2 

0750 

4.49 

— 

220.3 

.3564 

227.4 

0755 

3.70 

— 

85.9 

.1361 

86.8 

0800 

2.62 

32.4 

30.2 

.0543 

34.6 

0805 

2.54 

24.9 

23.2 

.0424 

27.0 

0810 

2.65 

31.8 

31.6 

.0593 

37.8 

0815 

2.66 

50.85 

63. 

.1155 

73.7 

0820 

2.83 

— 

87.4 

.1523 

97.1 

0825 

3.20 

92.3 

.1590 

101.4 
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18  FEB  81 0455 


Figure  18. 


Particle  Size  Distributions  for  Dense  Fog 
19  February  1981.  0455-0510  hrs 
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At  this  t'^me  a series  of  measurements  were  made  for  several  wavelengths. 
The  data  is  tabulated  in  Table  14.  Shortly  after  measuring  the  R16  line, 
the  signal  was  lost  and  was  not  detectable  during  the  next  40  minutes. 

The  s+cble  nature  o1  the  fog  between  0620  and  07i0  is  illustrated  in 
Figutc.  20-  This  is  reflected  by  the  good  correlation  between  o , a , 
and  a . Figure  21  presents  size  distributions  taken  ten  minutes  apart. 

As  car  be  seen  from  Table  13,  the  size  distributions  correspond  to  ex- 
tinction of  173.3,  270.9,  281  , and  220.3  d8/km..  The  highest  extinction 
computed  during  this  period  was  291  dB/km  and  occurred  at  0735  hours. 

Size  distributions  corresponding  to  extinction  data  taken  between  0805 
and  0820  hours  are  plotted  in  Figure  22.  The  fog  reached  its  n.inimum 
intep'ity  at  0805  hours  and  had  started  to  grow  in  intensity  at  0810 
hours.  The  part  of  the  distribution  with  radius  less  than  six  micro- 
meters is  almost  identical  to  the  same  part  of  the  distribution  obtained 
at  0805  hours.  The  small  incre  -'e  shown  in  the  larger  particles  between 
6 and  23  micrometers  is  "esp-  s'-hle  for  more  than  30%  rise  in  extinction. 

A similar  analogy  may  be  meas  for  the  distributions  shown  for  0815  and 
0820  'lOurs,  The  most  interesting  feature  of  Figure  22  is  that  the  mean 
and  maximum  radius  of  '.rr.'  distributions  shown  are  nearly  equal  in  size 
to  each  other,  yet  there  is  almost  a 4 to  1 difference  in  extinction 
produced.  Variation  in  extinction  with  changing  fog  is  documented  in 
Tdb.e  15.  Data  was  collected  with  the  R18  line  (10.26  ym).  A thirty- 
one-fold  decrease  in  transmission  is  observed  between  0808  and  0810  hours. 
Table  16  shows  transmission  and  extinction  data  for  several  CO^  laser 
lines  in  medium  density  fog.  Again,  it  is  seen  that  extinction  in  fog 
is  not  wavelength  dependent  even  for  medium  density  fogs.  The  larger 
increase  between  the  R18  and  the  R16  is  «.he  last  line  to  appear  in  the 
measurement  sequence,  and  by  then  the  fog  properties  have  slightly  changed 
along  the  optical  path. 

4.  RAIN 

Rain  c.vMnction  measurements  were  obtained  for  single  lines  as  a 
function  of  rain  rates  over  several  intervals  and  for  series  of  lines 
over  very  short  intervals.  Figures  23,  24,  and  25  are  plots  of  extinction 
versus  rain  rales  for  the  22pd,  28th,  and  29th  of  April  1981...  The  solid 
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TABLE  14 

EXTINCTION  OF  CO^  LASER  RADIATION  IN  HEAVY 

FOG  FOR  SEVERAL  WAVELENGTHS,  19  FEBRUARY  1931, 
0545-0546 

LINE  NO,  X (ya)  <r(dB/km) 


P28 

10.675 

$3.6 

P26 

10.653 

69.0 

P18 

10.571 

69.8 

P14 

10.532 

69.8 

R18 

10.260 

70.8 

R16 

10.275 

71.7 

Teaperature 

- 53.5*F 

Relative  Humidity 

Dew  Point  - 

53.5’F 

Pressure  ■■  990  ab 

1001 
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18  FEB  81  0805 


Figure  22.  Particle  Size  Distributions  for  Selected  Fogs, 
19  February  1981,  0805-0820  hrs 
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TABLE  15 


EXTINCTION  OF  10.26  ym  CO2  LASER  RADIATION 
WITH  TIME,  19  FEBRUARY  1981 


TIME 

T (10‘3) 

^ Oarh 

0'(dB/kin) 

0801 

9.44 

4.611 

20.0 

0802 

9.19 

4.69 

20.4 

0803 

6.00 

5.12 

22.2 

0804 

5.41 

5.22 

22.7 

0804 

3.24 

5.73 

24.9 

0805 

2.51 

5.49 

26.0 

0806 

5.08 

5.28 

22.9 

0807 

10.60 

4.54 

19.7 

0808 

20.60 

3.88 

16.9 

0809 

2.50 

5.99 

26.0 

0810 

.66 

7.32 

31.8 
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TABLE  16 

MEASURED  TRANSMISSION  AND  EXTINCTION  FOR 
SEVERAL  CO2  laser  LINES  IN  MEDIUM  DENSITY 

FOG,  19  FEBRUARY  1981,  0811-0812  HOURS 


LINE  NO. 

X (ym) 

T (10’^) 

^ (dB/km) 

P28 

10.675 

6.89 

31.6 

P26 

10.653 

6.68 

31.8 

P20 

10.591 

5.72 

32.4 

P18 

10.571 

5.46 

32.6 

P14 

10.532 

6.84 

31.6 

R18 

10.260 

5.95 

32.3 

R16 

10.275 

3.72 

34.3 

Temperature  - 52.2*F  Relative  Humidity  ■ 1007. 
Dew  Point  ■ 52. 2 "F  Pressure  " 990  mb 
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PRECIPITATION 


Figure  25.  CO.  Laser  Extinction  as  a Function  of  Precipitation 

L 

Rate  for  10.6  ym  Laser  Radiation,  29  April  1981 
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lines  are  a best  fit  model  obtainea  by  applying  a nonlinear  regression 
scheme  to  the  measured  data.  Data  obtained  for  April  22nd  contains 
extinction  values  for  rain  rates  up  to  58  mm/hr,  data  for  April  28th 

contains  extinction  values  for  data  up  to  A6  mm/hr,  and  data  for  April  29th 

contains  extinction  values  for  data  up  to  75  mm/hr.  The  total  number  of 

points  shown  can  be  misleading,  since  points  with  the  some  coordinates 
will  fall  on  top  of  each  other.  For  example,  in  the  data  shown  in 
Figure  25,  five  extinction  values  are  plotted  for  rain  rates  near  75  mm/hr, 
yet  only  two  data  points  are  discernible  at  this  value.  The  data  illustrated 
on  each  of  the  plots  was  accumulated  over  several  hours  of  measurement. 
Dynamic  situations  such  as  the  moving  in  and  moving  out  of  rain  showers 
were  observed  but  not  recorded.  Data  was  only  recorded  when  the  monitored 
signal  indicated  a semi-stable  optical  path  between  transmitter  and 
receiver. 

The  data  obtained  for  April  29th  contain  the  highest  rain  rates  and 
also  the  maximum  number  of  measured  points  (62  points).  A best  fit  curve 
to  the  measured  data  gave  the  relationship 

(T  - l.A33P®*^^^dB/kni  (141) 

where  P is  the  precipitation  rate  in  mm/hr.  The  data  for  22nd  and  28th 
April  were  combined  with  data  for  the  29th  and  the  result  is  plotted  in 
Figure  26.  Application  of  a nonlinear  regression  scheme  yielded  the 
equation 

(T  - 1.4P°*^dB/kni  (142) 

The  same  regression  scheme  was  applied  to  rain  extinction  data  published 
by  Chen  (Reference  54)  and  by  Rensch  and  Long  (Reference  55).,  The  re- 
sulting curves  are  plotted  alongside  the  curve  given  by  Equation  142. 

Chen's  calculations  were  obtained  from  data  derived  from  laws  and  Parsons 
rainfall  measurements.  The  rain  rates  were  derived  from  measurements  of 
drop  size  and  drop  terminal  velocity  (Reference  56).;  These  quantities 
do  not  lend  themselves  to  easy  measurements.  The  relationship  derived 
from  Chen's  results  is  given  by 

- 1.085P^'^^^dB/kro  (143) 
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while  that  derived  from  Rensch  and  Long's  results  is  given  by 

Cr^  - 1.84P®'^°^dB/km  (144) 

Equation  143  was  derived  from  data  with  rain  rates  up  to  50  mm/hr.  The 
extension  to  75  mm/hr  was  done  by  the  computer,  based  on  the  model  given 
by  Equation  143.  Data  from  Rensch  and  Long  was  obtained  from  direct 
precipitation  measurements.  Their  results  fo‘'low  nicely  those  of  Chen's 
up  to  rain  rates  of  25  mm/hr,  but  then  it  begins  to  drop  off.  Extinction 
values  as  a function  of  rain  rates  for  several  wavelengths  are  tabulated 
in  Table  18  for  a rain  rate  of  1,82  mm/hr  for  data  taken  on  April  22nd 
and  in  Table  19  for  rain  rates  of  1.8,  8.5,  and  12,0  mm/hr  for  oata  taken 
on  April  29th.  Whatever  differences  are  observed  are  due  to  the  differences 
in  molecular  contributions  and  to  changes  in  the  atmospheric  path  during 
the  time  of  measurement. 

5.  HAZE 

Due  to  the  uncertainty  in  the  aerosol  refractive  index,  extinction 
calculations  yield  very  uncertain  results.  Depending  on  the  relative 
humidity,  dry  aerosols  can  grow  up  to  ten  times  their  normal  size 
(Reference  57),  The  size  distribution  given  by  Junge's  model,  with  the 
calculated  growth  rate  in  humidity  is  shown  in  Figure  28.  Three  separate 
distributions  are  shown  for  three  different  type  atmospheres.  Curve  1 
is  for  a continental  atmosphere.  Curve  2 for  an  oceanic  atmosphere,  and 
Curve  3 is  for  the  average  sea  spray  (NaCl).  For  continental  aerosols, 
it  is  very  obvious  that  as  the  particle  grows  in  size,  its  refractive 
index  changes.  For  water  soluble  aerosols,  the  refractive  index  has  been 
measured  to  be  approximate  to  that  of  NaCl  (Reference  58),  As  the  parti- 
cle grows  with  humidity,  the  i^efractive  index  varies  in  values  to  that  in 
between  NaCl  and  H^O,  Therefore,  computations  of  extinction  using  Mie 
Scattering  theory  yie^d  very  rough  approximations  unless  the  exact  value 
of  the  aerosol  refractive  index  is  known.  This  is  illustrated  in  Tables 
19  and  20.;  Table  19  corresponds  to  the  aerosol  distribution  shown  in 
Figure  29  and  Table  20  corresponds  to  Ihe  aerosol  distribution  shown  in 
Figure  30.  The  aerosol  extinction  was  computed  usirig  a refractive  index 
approaching  that  of  NaCl  (n  = 1 .45  - .154  j)  duu  .hat  of  H^O  (n  = 1.185  - 
.0662  j)„  Comparing  the  extinction  due  to  water  soluble  aerosols  (o^) 
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TABLE  17 

EXriNCTION  OF  CO2  LASER  LINES  FOR  A CONSTANT 
RAIN  RATE  OF  1.82mm/HR,  22  APRIL,  1935  HOURS 


LINE  NO. 

X (ym) 

^ (dB/Vm) 

P28 

10.675 

2.78 

n 

P22 

10.611 

2.72 

P20 

10.591 

2,67 

P18 

i 

10.571 

2.69 

P14 

10.532 

2.68 

, 

R16 

10.275 

2.82 

R18 

10.260 

2.96 

Temperature 

- SOT 

Relative  Humidity  - 99% 

Dew  Point  ■ 

59. 5T 

Pressure  ~ 990  mb 
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TABLE  18 

EXTINCTION  OF  SEVERAL  CO2  LASER  LINES  FOR 

CONSTANT  RAIN  RATES  OF  1.8,  8.5,  and  12  mm/HR, 
20  APRIL  1981 


LINE 

NO 

X 

ym 

EXTINCTION 

P - 1.82  P - 8.5 

P - 12 

P28 

10.675 

2.78 

5.65 

5.86 

P22 

10.611 

2.74 

5.47 

6.21 

P20 

10.591 

2.69 

5.69 

6.18 

P18 

10.571 

2.61 

5.91 

6.12 

P14 

10.532 

2.65 

5.73 

5.99 

R16 

10.275 

— 

5.73 

5.73 

R18 

10.260 

3.04 

6.30 

I 
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TABLE  19 

COMPARISON  OF  AEROSOL  EXTINCTION  HAVING  AN 
ASSUMED  COMPLEX  INDEX  OF  H2O  AND  AN  ASSUMED 

COMPLEX  INDEX  OF  NaCi  (DATA  TAKEN  27  JULY 
1981) 


ym 

km"l 

<r. 

Vcm”! 

to-1 

AT 

km“l 

.46 

.088 

.242 

.585 

.533 

.052 

.49 

.090 

.247 

.588 

.533 

.058 

.46 

.115 

.316 

.588 

.533 

.058 

Temperature 

- 73  *F 

Relative  Humidity 

- lOOX 

Dew  Point  ■ 

69.5*F 

Absolute  Humidity 

- 17.5 
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Figure  27.  Best  Fit  Models  for  Data  on  Figure  26,  Laws  and 
Parse  ns  as  Reported  by  Chen,  and  Measurements 
Conducted  by  Rensch 
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10*  10*  10^  10®  iir  M 


RADIUS 


Figure  28.  Model  Size  Distribution  for  Atmospheric  Aerosols, 
Adapted  from  Reference  27,  p.  124 

The  hatched  area  between  Curves  2 and  3 represents 
the  non- sea  spray  component  over  the  ocean.  The 
straight  lines  indicate  the  shift  of  the  continental 
log  radius  number  distributions  curve  due  to  particle 
growth  by  humidity  for  mixed  particles  of  about  20% 
soluble  material . 
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TABLE  20 

COMPARISON  OF  AEROSOL  EXTINCTION  HAVING  AN 
ASSUMED  COMPLEX  INDEX  OF  H^O  AND  AN  ASSUMED 

COMPLEX  INDEX  OF  NaCl  (DATA  TAKEN  5 AUGUST 
1981) 


% 

ViQ 

km’^ 

lan-l 

kai"l 

CTp 

km-1 

ACT 

laa“l 

.43 

.049 

.135 

.536 

.498 

.038 

.42 

.042 

.132 

.535 

.498 

.037 

.42 

.032 

.096 

.531 

.498 

.033 

Temperature 

- 68.4'F 

Relative  Hxnnidlty 

- 88 

Dew  Point  ■ 

68.4'F 

Absolute  Humidity 

- 18 

gr/m 
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to  that  of  water  aerosols  (o  ),  it  is  seen  that  a difference  in  extinction 

w 

of  2.7  to  1 , or  greater,  exists  in  all  cases.  The  aerosol  extinction 
values  were  computed  for  a wavelength  of  10.6  microns.  The  difference 
between  the  molecular  extinction  and  the  measured  extinction  is  given 
by  La. 

Two  factors  must  be  remembered  in  this  comparison.  The  most  important 
one  is  that  a very  few  number  of  large  particles  (4  or  5)  can  unduly 
influence  the  aerosol  extinction  calculated.  This  is  clearly  shown  in 
Table  21.  If  the  number  of  particles  measured  with  radius  3.373  ym  would 
have  been  five  instead  of  one,  the  extinction  for  that  size  radius  would 
be  0.009  instead  of  0.0018.  The  missing  of  only  two  particles  of  radius 
3.62  ym  results  in  the  omission  of  0.045  from  the  aerosol  extinction 
computed.  As  can  be  seen  from  Table  21,  particles  with  size  radius  less 
than  0.35  ym  do  not  contribute  significantly  to  the  aerosol  extinction. 
Although  the  lower  number  density  has  a larger  number  of  particles,  the 
significant  difference  to  the  contribution  of  extinction  lies  almost 
entirely  to  the  small  Increase  on  the  density  of  particles  with  mean 
radii  of  0.612  and  0.866  microns.  The  second  factor  Is  due  to  the  error 
uncertainty  of  the  particle  measuring  system  in  measuring  sizes  at  the 
low  range  of  the  Instrument.  A third  possible  source  of  difference  could 
be  in  the  computed  value  of  given  by  the  FASCODE  model.  As  mentioned 
earlier  in  Section  III,  the  molecular  extinction  on  humid,  summer  days 
is  due  mostly  to  the  contribution  of  the  water  vapor  continuum.  The 
exactness  and  source  of  the  contribution  of  the  water  vapor  continuum 
is  still  being  debated. 

Measured  transmission  values  are  compared  with  calculated  values 
for  several  wavelengths  In  Tables  22  and  23.  Aerosol  size  distributions 
corresponding  co  extinction  data  of  Tables  24  and  25  are  plotted  in 
Figures  29  and  30.  The  difference  in  transmission  for  the  P^q  line  Is 
0.039  for  the  0940.  July  27th  measurement.  However,  the  loss  In  trans- 
mission due  to  aerosol  extinction  (a  = .09)  is  approximately  0.10  and 

d 

this  leaves  a factor  of  0.06  unaccounted  for.  A plot  of  the  FASCODE  out- 
put corresponding  to  the  data  tabulated  in  Table  23  is  shown  In  Figure  31. 
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TABLE  21 

COMPUTATION  OF  EXTINCTION  FOR  PARTICLES  OF 
RADIUS  r(I)  AND  DENSITY  dN(r)/dr  (5  AUGUST 
1981,  0800  HOURS 


r(I) 

N(I) 

dN(r)/dr 

^(i) 

.353 

8132 

71.900 

.0010 

.612 

2409 

21.300 

.0249 

.866 

363 

3.210 

.0093 

1.118 

13 

.115 

.0005 

1.369 

7 

.062 

.0005 

1.620 

6 

.053 

.0008 

1.871 

5 

.044 

.0011 

2.121 

4 

.035 

.0014 

2.372 

4 

.035 

.0021 

2.622 

3 

.027 

.0022 

2.872 

2 

.018 

.0020 

3.122 

1 

.009 

.0013 

3.373 

1 

.009 

.0018 
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TABLE  22 

MEASURED  AND  COMPUTED  TRANSMISSION  FOR 
SEVERAL  CO2  LASER  LINES  IN  HOT.  HUMID, 

HAZY  WEATHER,  27  JULY  1981,  0940  HOURS 


LINE 

NO. 

X 

pm 

Tjn 

(measured) 

Tc 

(computed) 

AT 

P28 

10.675 

.571 

.607 

.036 

P26 

10.653 

.566 

.602 

.036 

P22 

10.611 

.556 

.590 

.034 

P20 

10.591 

.548 

.587 

.039 

P18 

10.571 

.554 

.590 

.036 

R16 

10.275 

.561 

.595 

.034 

R18 

10.260 

.573 

.608 

.035 

Temperature  - 20.21'C 
Dew  Point  - 20.21*C 


Relative  Humidity 
Absolute  Humidity 


lOOX 

17.5  3 
gr/m-^ 
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TABLE  23 

MEASURED  AND  COMPUTED  TRANSMISSION  FOR 
SEVERAL  CO^  LASER  LINES  IN  HOT,  HUMID, 

HAZY  WEATHER,  5 AUGUST  1981,  0825  HOURS 


Line 

No 

B 

(measured) 

Tc 

(confuted) 

AT 

P28 

10.675 

.602 

.622 

.020 

P26 

10.653 

.596 

.618 

.022 

P22 

10.611 

.590 

.613 

.023 

P20 

10.591 

.584 

.608 

.024 

P18 

10.571 

.585 

.610 

.025 

R16 

10.275 

.593 

.613 

.020 

R18 

10.260 

.604 

.623 

.019 

Temperature  - 22.75*C  Relative  Humidity 
Dew  Point  * 20. 8“ C Absolute  Humidity 


88.8% 

18.8  , 
gr/m'5 
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FOR  SEVERAL  CO^  LINES. 

THE  INABILITY  TO  OBTAIN 
SHORT  TIME  REQUIRED 

T2  T3 


.051 

.051 

.225 

.049 

.186 

.049 

.189 

.047 

.187 

.048 

.190 

.044 

.173 

I 

! 
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6.  SNOW 

No  reliable  means  exist  at  our  facility  for  categorizing  snow. 

Although  visibility  is  not  a good  tool  for  this  purpose,  it  is  the  only 
one  at  our  disposal.  Three  sequences  of  transmission  measurements  through 
snow  are  shown  in  Table  24.  The  three  transmissions,  T^ , T^,  and  T^ 
correspond  to  approximate  visibilities  of  1.5,  1.2,  and  1.7  km,  respectively. 
All  three  were  taken  under  what  was  considered  light  snowfall  condition. 
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SECTION  VII 
CONCLUSIONS 

Results  show  that  extincitlon  of  CO,  laser  radiation  in  dense  fog 
may  be  determined  to  a good  approximation  by  the  relation  a » 147W,  where 
W is  the  liquid  water  content  of  the  fog  distribution.  With  few  exceptions, 
results  obtained  by  use  of  this  approximation  give  results  which  are 
within  7%  of  results  obtained  by  use  of  Mie  Scattering  theory.  The 
approximation  seems  to  be  good  for  extinction  values  greater  than  9 dB/km. 

Extinction  by  dense  fog  can  easily  exceed  the  dynamic  range  of  the 
measuring  Instrument.  The  probability  of  this  occurring  is  better  than 
0.5.  Use  of  higher  sensitivity  detectors  does  not  appreciably  change 
this  number  when  dealing  with  extinction  in  excess  of  70  dB/km. 

Extinction  by  rain  with  precipitation  rates  up  to  75  mm/hr  was 
measured.  Appl ication  of  an  empirically  derived  model  gives  good  results 
with  data  obtained  by  other  workers.  The  derived  model  can  be  used  to 
obtain  extinction  from  knowledge  of  the  precipitation  rate  without  need 
for  measuring  the  rain  drop  size  and  terminal  velocity.  Measured  extinc- 
tion in  heavy  rainfall  was  found  to  be  an  order  of  magnitude  less  than 
measured  extinction  in  heavy  fog. 

Extinction  by  summer  aerosols  is  difficult  to  calculate  due  to  the 
uncertainty  on  the  composition  of  the  aerosols.  An  upper  and  a lower 
bound  on  the  aerosol  extinction  has  been  calculated  using  the  complex 
index  of  NaCl  and  H^O,  respectively. 

Comparison  of  measured  distributions  to  those  given  by  models  reveals 
that  these  models  are  only  applicable  for  those  cases  where  smooth  con- 
tinuous distributions  are  measured.  Unfortunately,  this  type  of  distri- 
bution was  seldom  observed,  and  in  most  instances,  smoothing  approximations 
can  only  be  applied  to  less  than  25*  of  the  measured  distributions.  No 
smooth  haze  distributions  were  observed.  This  indicates  that  models  for 
size  distributions  can  give  serious  errors  when  used  in  computing  aerosol 
extinction  and  should  be  used  cautiously. 
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Extinction  measurements  for  CO^  laser  propagation  under  adverse 
weather  clearly  points  out  that  extinction  is  at  best  a very  weak  function 
of  wavelength.  This  can  be  seen  by  looking  at  the  extinction  cross  section 
Q(x)  that  depends  on  the  ratio  of  X/r  of  the  Mie  size  parameter  x = 2Tr\/r. 

For  CO2  laser  wavelengths,  the  departure  from  10.6  ym  is  not  that  significant 
in  terms  of  wavelength  (±0.4  ym),  but  the  true  wavelength  dependence  of 
the  extinction  coefficient  has  to  be  examined  in  terms  of  the  complex  re- 
fractive* index  of  water.  For  clouds  with  a mode  radius  of  4.0  ym, 
Oeirmendjian  gives  the  extinction  at  10.0  ym  to  be  11.18  km"^  and  10.10  km'^ 
at  11.5  ym  (see  p.  189  of  Reference  34).  This  is  approximately  a ±5* 
deviation  from  the  10.6  ym  CO^  wavelength  of  interest. 

A final  conclusion  is  that  a CO2  laser  radar  with  a 3 dB/km  performance 
has  ve.y  limited  utility  in  fog  (except  for  light  winter  fogs),  may  be 
used  in  rain  with  precipitation  rates  of  25  mm/hr  at  ranges  of  1 kilometer 
or  less,  and  has  little  utility  in  snow  except  for  very  light  snowfall 
conditions . 
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APPENDIX  A 

OORRECnON  FACTORS  TO  THE  VAN  de  Hulst 
APPROXIMAnONS  TO  THE  MIE  EXTINCnON  POIMJLA 

The  Deitmendjlan  oorrectlan  factors  to  the  Van  de 
Hulst  appiaximtlons  are  given  by  (see  Reference  21,  page 
30). 

(n<,  - 1)2  0.2p  - n-i  + 1 

D,  - (S  + 11  + 1 , p - 5(n,  - 1) 

^ 1.632n,  (n,  - 1)S  ^ 


CA-1) 


tu  - 1 
• _J (S  + l)p  , 


8.16n, 


5(ni  - 1)  ^ p ^ 


4.08 


1 <1*  3tan0 


(A-2) 


0*1-  1)(S+  1) 

2ni<l  + 3tanS)  1+3  tang  1 + tmS 


4.08  ^ 4.08 

S < 


P - 


(A-3) 


2.04(v  - 1)<S  + 1) 


4.08 


P > 


1 + toiB 


CA-4) 


\here 


S » 1 + 4tan6  + 3 tan^B 


(A-5) 


niMm 
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p “ 4Trr/X(n]^  - 1)  (4-6) 

B “ t:an"^(nj^/(n2  “ 1)1  (4-7) 

X « 27Tr/  - 1)  (4-8) 

For  = 1.185,  = .0662  and  X = 10.6  microns,  Equations  A-5  through 

A-8  reduce  to 

5 - 2.816 

p - 0.2193r 

6 - 0.358 

X - 0.0785r 


v;here  the  radius  r is  expressed  in  microns.  Using  the  above 
values,  it  is  found  that  the  correction  factors  and  their 
ranges  are, 


0.0675  + 1.92(.0439r  - 

0.185)  r - 4.21 

(A-9) 

o.ieir 

4.21  < r ^ 8.97 

(A-10) 

.1434 

8.97  < r ^ D.7 

(A^ll) 

.5158/r 

r > 13.7 

(A-12) 

Ihe  extinction  value  is  obtained  from 
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Use  of  Equations  A-10  through  A-12  in  conjunction  with  Equation  85  gives 
extinction  values  which  are  within  6%  of  exact  values  obtained  from  Mie 
formula  calculations.  For  r = 0.577  microns,  use  of  gives  an  error 
of  Z7%,  and  at  the  upper  end  of  its  useful  approximation  it  gives  a 6.U 
error  for  r = 3.51  microns. 

(^  - Qg(l  + D)  (A-13) 

In  order  to  reduce  the  error,  the  Ueirmendjian  correction  factors 
were  modified  to  give  better  approximations  when  used  with  the  Van  de 
Hulst  equation  for  extinction.  To  achieve  this,  it  was  necessary  to 
break  down  the  approximation  factors  in  the  range  0 < r £ 4.0  into  two 
different  factors  with  ranges  0 < r < 1 and  1 £ r < 4.0.  The  correction 
factors  and  their  useful  range  are  given  below,  where  every  attempt  was 
made  to  retain  the  original  form  of  the  Deirmendjian  factors. 


(ni  - 1)2  (S  + 4.6) 


1.44n, 


(n^  - 1)  - pr" 
5(n^  - 1)S 


r < I 


CA-IA) 


Oai  - 1)^  (S  + Dr' 
3,673n, 


- 1) 


(In  r)' 


(n^-l) (s+l)pr‘ 


165n, 


I - r < 4.  0 ^-15) 


4.0  - r - 7.0 


(A-16) 


1.794(nj-l) 

(l+3tanB) 


7.0  i r 1 13.7 


(A- 17) 


1.89(nj-l)(S+l) 
S tn  r 


r i 13.7 


(A- 18) 


127 


AFWAL-TR-81-1280 


Use  of  these  factors  in  conjunction  with  Equation  A-13,  give 
values  which  are  within  M of  values  obtained  from  exact  Mie 


extinction 

calculations. 
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APPENDIX  B 

CONVERSION  OF  EXTINCTION  UNITS 
From  km  ^ to  dB/km 


The  extinction  coefficient  o can  also  be  expressed  in  decibels  (dB) 

A 

by  use  of  the  standard  definition 

dB  - 10  log  (Iq/I) 

- 10  log  (1/T) 


- 10  log  0 (B-1) 

where  0 is  the  opacity  of  the  medium.  The  optical  density  of  the  medium 
is  related  to  opacity  by  the  relation 

D - log  (1/T)  (B-2) 

or  .n 

T - 10  ^ (B-3) 


Substituting  Equation  B-2  in  Equation  B-1  one  obtains  the  relation 


dB  - lOD 

or 

D - .IdB  (B-4) 

Using  a unit  path  length  { tinctlon  coefficient  has  units  of  km’^),  the 
atmospheric  transmittance  given  by 

(Tx 

T - c*'  ^ (B-5) 

Comparing  Equation  B-3  and  Equation  B-5  ore  sees  that 


e 


(B-6) 
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Taking  logs  of  both  sides  of  Equation  B-6  one  obtains, 

* D In  10 


2.302585  D 


Substitution  of  Equation  B-4  in  Equation  B-7  yields 


(f  « 0.230585  db 

A 


(B-7) 


(B-8) 


The  attenuation  loss  per 


kilometer  is  easily  obtained  from  Equation  B-8 


f(dB/km)  - 4.3367936  C 
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APPENDIX  C 
TEST  SITE 

Because  of  safety  considerations,  the  propagation  path  had  to  be 
located  in  a controlled  area  with  a requirement  for  beam  termination  on 
the  ground  in  case  the  laser  beam  missed  the  receiving  van.  The  site 
chosen  was  adjacent  to  an  old  runway  on  Area  B,  Wright-Patterson  AFB 
(Figure  C-1).  The  receiver  was  positioned  of^  the  runway  so  as  to  eliminate 
the  thermal  heat  emanated  from  the  cement  structure.  Because  the  wind 
normally  blows  from  the  southeast,  the  heat  from  the  runway  would  often 
be  carried  into  the  optical  path  creating  extreme  turbulence  fluctuations 
which  had  to  be  dealt  with  (see  Section  V).  Because  of  the  isolated  site, 
diesel  power  generators  were  used.  These  generators  were  very  unreliable, 
and  at  times  impossible  to  start  when  the  temperature  fell  below  the 
freezing  point.  Since  the  generators  could  not  be  left  running  continu- 
ously, the  heaters  in  the  receiver  van  were  normally  off.  This  resulted 
in  condensation  forming  on  the  24-inch  mirror  whenever  the  dew  point  was 
high  and  the  mirror  was  colder  than  the  surrounding  air.  Whenever  this 
happened,  it  was  necessary  to  wait  a few  days  until  the  mirror  temperature 
rose  higher  than  the  surrounding  air  temperature  before  cleaning  of  the 
mirror  was  possible.  Eventually,  gas  heaters  were  installed  and  the 
problem  resolved. 
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